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Abstract 

We derive factorization theorems for Aqcd/^ti;, power corrections to inclusive B-decays in the 
endpoint region, where ~ mi, Aqcd . In i? — > Xutv our results are for the full triply differential 
rate. A complete enumeration of Aqqd/ mi, corrections is given. We point out the presence of new 
AqcD/mfc-suppressed shape functions, which arise at tree level with a 47r-enhanced coefficient, and 
show that these previously neglected terms induce an additional significant uncertainty for current 
inclusive methods of measuring \Vub\ that depend on the endpoint region of phase space. 
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I. INTRODUCTION 



The study of -B-decays enables precision measurements of Standard Model parameters 
as well as searches for new physics. For inclusive measurements, the large mass of the b- 
quark makes a rigorous theoretical treatment in QCD possible. A well-known example is 
the semileptonic decay B X^iiy, which allows measurements of \Vcb\, fnh and rric through 
moments of the decay spectra. Analyses of this type are performed at BaBar Belle 
and CLEO 0], where simultaneously fitting_a few fundamental hadronic parameters keeps 
the theoretical uncertainties under control p. From measurements of -B — >■ we can 



also determine \Vub\ |y,|7|, |8|. Finally, in channels involving flavor-changing neutral currents, 
such as B Xg'-f and B — > Xsi~^i~ , we have the luxury of searching for new physics with 
accurate Standard Model calculations of the expected rates 

In inclusive S-decays there is often a trade-off between theory and experiment, because 
cuts are necessary experimentally, but these less inclusive spectra make the theory more 
complicated. For sufficiently inclusive measurements we can use quark-hadron duality, to- 
gether with a local operator product expansion (OPE) in AqcB/mb Hil, 13 H H H E3 ■ 



In these decays, non-perturbative matrix elements are defined with the help of the Heavy 
Quark Effective Theory (HQET) [13,0, and perturbative corrections are straightforward 
to incorporate, although often still difficult to compute. For B — > JC£u decays, this formal- 
ism has been worked out to order l/mf in the power expansion |l9|,|20| and the perturbative 



corrections are in most cases known to a^/^o 122, |23|, |2J| and in a few cases at or- 
der a?_25j_2^|. A simultaneous fit with the hadronic parameters appearing up to order 
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1/ml |27il28l| has provided us with a ~ 2% determination of \Vcb\, making it one of the best 
determined elements of the CKM matrix. 

In B ^ Xulv decays the situation is more complicated because of the large b ^ c 
background. Since a model-independent analysis must avoid extrapolating the experimental 
data, the computation of decay rates in a restricted region of phase space becomes important 
(see e.g. 0, E3])- Phase-space cuts are also important for B — > Xg'y to eliminate softer 
photons, and for B Xsi^i^ to deal with the intermediate J/"^ and resonances. In 
many cases, these restrictions on the phase space can put us in a situation where the local 
operator product expansion is no longer applicable, since we are limited to a region where 



the allowed masses of the hadronic states satisfy m'^ -C ml. To distinguish different cases 



we use the terminology 

~ m^, totally inclusive (local OPE, HQET) 



Arrix ~ "^feAqcD, endpoint region (Factorization, SCET) , (1) 

2 

QCD' 



Am^ ~ ^QCD' resonance region (exclusive methods) 



where Am^ denotes the region of extending out from iTi'j^^.^, and the applicable theo- 
retical methods are shown in brackets. 

When 171^ ~ mftAqcD, the outgoing hadronic states become jet-like and we enter the 
so-called endpoint or shape function region [i3,E3, 0, which will be the main focus of this 
paper. For h ^ u decays this region is important, because of cuts on or m^, which are 
used to eliminate b —>■ c events. In this region both the perturbative expansion and the 
power expansion become more complicated. In particular, there is the usual perturbative 
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expansion at the scale /i^ ~ ml, as well as a second perturbative expansion at the smaller 
scale fi^ ~ mj^. The rates also exhibit double Sudakov logarithms. In addition, instead 
of depending on non-perturbative parameters (Ai, A2, . • .) that are matrix elements of local 
operators, the decay rates depend on non-perturbative shape functions. We shall refer to 
the expansion parameter for this region as ~ m\/m\ ~ Aqcd/'"^?) to distinguish it from 
the l/m^ expansion for the local OPE. In the endpoint region, the standard OPE no longer 
completely justifies the separation of short- and long-distance contributions. Instead, we 
must consider a more involved derivation of a QCD factorization theorem, as is the case in 
processes such as Drell-Yan and DIS as x — > 1 [j^. l35|. 

For B —>■ Xsj and B —>■ Xuii' at leading order fLO) in A, the factorization theorem for 
the endpoint decay rates was determined in Ref. It separates QCD contributions that 
are hard (H), coUinear {J^^^^) and soft (Z^^-*), so that, schematically, a differential decay rate 
takes the form 

dr = Hx j^(o) ® /(°) , (2) 

where x is normal multiplication and is a one-parameter convolution. Here the hard 
contributions are perturbative at the scale fj? ~ ml, the coUinear contributions in J'^^^ are 
associated with the inclusive X jet and are treated perturbatively at the scale /i^ ~ m^AqcY), 
and the soft contributions are factored into a forward S-meson matrix element giving the 
non-perturbative shape function /^"^ js^- In Ref. js^] this factorization theorem was red- 
erived using the Soft-Collinear Effective Theory (SCET) El H 0. The attraction 



of the effective-theory method is that it provides a formalism for extending the derivation 
of factorization theorems beyond LO in the power expansion. The main goal of this paper 
is to derive a factorization theorem for B — »• Xgj and B —>■ X^iu at subleading order, 
i.e. C(A^), using SCET. This factorization theorem allows us to separate perturbative and 
non-perturbative corrections to all orders in a^. 

One method for studying the endpoint region is to start with the local OPE and sum up 
the infinite series of the operators that are most singular as we approach the ~ ^b-^QCB 
region. This technique was used in Refs. j3l[ 3^ 3^, 41,|4^|43[, and provides a method of 



defining the non-perturbative functions. At LO the result is the shape function^ 



S^^\t) {B^\h^{m-D)''K\B^) = -{B^\KS{t -in-D)hy\By) 

(3) 



A:=0 



The result is simply the matrix element of a non-local HQET operator, where the states \By) 
and heavy b-quark fields are defined in HQET, and is a light-like vector along the axis 
of the jet. This approach allows direct contact with the extensive calculations made with the 
local OPE, which give terms in the power series. There are several reasons for considering 
an approach where f{i^) is obtained without a summation. In particular, it is difficult to 



We arbitrarily use the term "shape function" for f{£'^) rather than "distribution function". Sometimes in 
the hterature the term "shape function" is reserved for the distribution that enters dr(B X^^v) / dEi, 
which is an integral over /. 
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go beyond tree level with the summation approach. Also, owing to the presence of a kinked 
Wilson line ji^], the renormalization of the local operators in the sum and final delta-function 
operator are not identical |45, (see also [l^), essentially because the moment integrals 
introduce additional UV divergences.^ For this reason the expansion in Eq. Q should be 
considered to be formal, and care must be taken in drawing conclusions from operators in the 
expanded version, such as the fact that they have trivial dependence on n^. Care must also 



be taken in calculating the hard factor H , as pointed out recently |45l. I46| . since the one- loop 
matrix element for f^^\^'^) has finite pieces in pure dimensional regularization. This implies 
that the quark-level QCD computation does not directly give the hard contribution, unlike 
factorization theorems involving massless quarks such as in DIS. The matching calculation 
in SCET handles this in a simple manner because matrix elements of the effective-theory 
graphs are necessarily subtracted from the full-theory graphs in order to compute H . 

For any precision calculation, perturbative corrections play an important role, and both 
the resummation of large logarithms and fixed-order calculations need to be considered. 
The position-space version of E g. Q has a kinked Wilson line along v-n-v, which leads to 
double Sudakov logarithms |36l. 144^. These occur between both the ml mbAqcD and 
the mfeAqcD ^qcd scales and can be summed using renormalization-group techniques. 
In moment space the leading and next-to-leading anomalous dimensions can be found in 
Ref. H. For phenomenological purposes, formulae for the differential rates with resummed 



logarithms are of practical importance and were obtained using inverse Mellin transforma- 



tions in Refs. |42l, |4Sl, |49l, l5(lj. These resummations have also been considered in SCET, 
both in moment space f38« . .4,^] and for the differential rates [i^ . Finite-order perturbative 
corrections are currently known to order for the H and functions 45. iol . 



Since in the endpoint region Aqcd/<5 ~ 1/5 to 1/10, it is important to consider the effect 
of power corrections for precision phenomenology. By matching from QCD at tree level, 
contributions of subleading NLO shape functions have been derived for B Xg'j and 
for B — > X^iu Is^lssl. followed by further analysis in Refs. jHi, 55, 5^. A single NNLO con- 
tribution has also been considered, cor resp onding to the "annihilation" contribution, which 
is phase space enhanced by IQn These power corrections provide the dominant 

source of theoretical uncertainty in current measurements of \Vub\ and are the focus of our 
paper, so we discuss them in more detail. 

To build some intuition, it is useful to contrast the power expansion in the endpoint 
region with the expansion for the local OPE. For the local OPE, all contributions can be 
assigned a power of (AQCD/''^fe)'^~^- The power of Aqcd and thus k is simply determined by 
the dimension of the operator, and the —3 accounts for the dimension of the HQET states. 



^ From the point of view of effective field theory this makes sense, since the summation in Eq. attempts 
to connect one effective theory (HQET) to a region with a different expansion parameter that is described 
by a different effective theory (SCET). Generically the renormalization in two EFTs is not interconnected. 
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FIG. 1: Comparison of the ratio of annihilation contributions to the lowest-order result. In the 
total decay rate, b) is ~ 167r^(A^/m^)Ai? ~ 0.02, while c) is ~ 47ras(m6)(A^/m^) ~ 0.003 when 
compared to a). In the endpoint region, b) is ~ 167r^(A^/m^)A5 ~ 0.16, a large correction, while 
c) becomes ~ 47ras(^j)(A/m5)e' ~ 0.6e', a potentially large correction. 

{By \ ■ ■ ■\By). For example, the set of local operators up to dimension 6 is 

0^ = KK, 0''^ = hSDTfK, 0''' = gKa^pG''^K, (4) 

= (7^.7"gL) {uiM , o'"' = {K qi) (jiL K) , 

O^^ = {KT'^rqL) {IlT^iM , O^^f = {KT'^qL) (jlL T"" K) , 

where dimensions are shown as superscripts, a superscript/subscript T means transverse 
to the HQET velocity parameter v^, and an L means left-handed.^ Dimension-4 oper- 
ators are absent so there are no l/m^ corrections, except the trivial ones that may be 
induced by switching to hadronic variables. For dimension-5 and 6 operators there are two 
naming conventions in common use. For l^B^\{0^°' ,0^^ ,0^°" ^0'^^}\B^) , the parameters are 
{Ai, A2, pi, P2} or {/i^, pI), p\s}- These operators are generated by connected graphs from 
the time-ordered product of two currents, as in Fig. QJi. On the other hand, the four-quark 
operators O^'^'^'^ give parameters f^Bi 2 and are disconnected (or rather connected by leptons 
or photons only), as shown in Fig. and thus exhibit a phase-space enhancement relative 
to Fig. The simplest way to see this is to note that for the total rate to B —>■ X^ii^, we 
would cut a two-loop graph for Fig. [T^, while Fig. [T|d would be at one-loop level (the i-~i> 
loop). For later convenience, we also consider the perturbative correction to the four-quark 
operators shown in Fig. [TJ:;, which is suppressed by as/(47r) relative to Fig. [TJd, and gives the 
operators O^^'^-^. In the total decay rate, if we normalize so that Fig. ^ ~ 1 then 

Fig. lb ~ IGtt^ ^ A5 ~ 0.02 , Fig. Ic ~ 4™ Jm;,) ~ O.OOSe . (5) 

Here AB = B2 — Bi 0.1 accounts for the fact that the matrix elements of the operators 
generated by Fig. ^ vanish in the factorization approximation. The factor of e accounts for 
any dynamical suppression of Fig. The definitions of B12 are 

{B,\ [K^^qL] [qilrK] \B,.) = [(fil - B2)gar + (4^2 - Bi)v^Vr ] . (6) 



We write in terms of HQET fields, although strictly speaking at lowest order this is not necessary. 
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Without the AB suppression factor, Fig. [TJd would dominate over other I /ml operators 
rather than just competing with them. The 0{as) corrections to annihilation are still a 
small contribution in the local OPE, for any e < 1. In particular possible enhancements of 
these contributions have been shown to cancel for the total 6 — >■ u decay rate js^. 

In the endpoint region there are extra enhancement factors and the dimensions of the 
operators no longer determine the size of their contributions. The fact that annihilation 
effects are larger in the endpoint was first pointed out in Ref. ji^. The power counting 
in SCET organizes these contributions in a systematic fashion and allows us to be more 
quantitative about how large these contributions are. Since some background material is 
required, we postpone this power counting until Sec. IIIII The derivation given here is more 
heuristic, but leads to the same results. For Fig. la the intermediate quark propagator 
becomes coUinear, giving an m^/A enhancement. This explains why a larger portion of the 
decay rate is concentrated in the endpoint region. For Fig. lb there is no quark-propagator 
enhancement but also no reduction from the phase space. A numerical estimate for this 
contribution was made in Ref. |57[. Finally, for Fig. Ic in the endpoint region there can be 
three coUinear propagators, giving a large mf/A^ enhancement to this diagram. In Sec. IVII 
we show that this graph contains the maximum possible enhancement. In summary, if we 
consider the rate integrated only over the endpoint region then Fig. la ~ 1 and 

A2 a 
Fig. lb ~ levr^ ^AB ~ 0.2 , Fig. Ic ~ 47rail.4 GeV) — ~ 0.6 e' . (7) 

nib 

The non-perturbative function that gives e' differs from the local operators that give e in 
Eq. (jSJ. Since e' ~ 0.3 is possible^, we conclude that the contribution from Fig. Ic may 
actually give one of the largest uncertainties in extracting Vub with methods such as Ei or 
m\ cuts that depend on the endpoint region. It has not been considered in recent error 
estimates in the literature. The main phenomenological outcome of our analysis is a proper 
consideration of this term for endpoint spectra. 

Theoretically, the main result of our analysis is a complete theoretical description for the 
NLO term, F^^), in the power expansion of decay spectra in the endpoint region. 



dT 
dZi 



+ -— + .... (8) 



endpoint 



Here Zi denotes a generic choice of the possible spectrum variables, 
{P^ , P~ , E^,q'^, SH,mb, . . .}. At NLO we use SCET to determine the contributions 
to the spectra. These contributions are tabulated in the body of the paper, but the generic 
structure of a term in (l/To)dT^'^^ /dZi is 

[dzn']H^^^\z^,,,mb, Zi) [ J^^'\zr.,, P±) , (9) 



^ The only rigorous scaling argument for a dynamical suppression that we are aware of is the large- iVc limit, 
where AB ~ l/Nc- For Fig^ the qq contraction gives a matrix element that is leading order in Nc- If 
these contributions are small, then a rough estimate is e' ^ l/Nc ~ 0.3. 
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where the number of convolution parameters varies from n = 1 to n = 3 and n' = 1 or 2, and 
for n = 2 [dk'^] = dkfdk2 etc. The dependence on the Zn' parameters appears only in jet 
functions that vanish at tree level. In Eq. (jU)) the (ji), (j2), (js) powers indicate whether the 
power suppression occurs in the hard, jet or soft regions respectively. The power corrections 
start at O(A^), which is ~ l/mb, and so ji + j2 + js = 2. Here ji,3 > while j2 can be 
negative. Phase-space and kinematic corrections give an 

if (2) with the same jet and shape 
functions as at leading order. Other more dynamic power corrections involve new hard H^^^ 
functions, and obtain their power suppression from the product of jet and soft factors. We 
show that the operators at NLO allow —4 < 72 < 2 and < js < 6. The largest jet 
function (^2 = —4) occurs for exactly the endpoint contribution generated by the four-quark 
operators (j3 = 6) from Fig. Ic. 

Our analysis can be compared with the closely related physical problem of deep inelastic 
scattering with ^ A^, in the limit where Bjorken x ~ 1 — A/Q. With no parametric 
sca ling for x, the power corrections in DIS at twist 4 were computed in Refs. |59[l6d 61, 6^, 



e3|, |6^. As a; — *■ 1 the relative importance of the power-suppressed operators changes and 
the importance of contributions from four-quark operators has been discussed in Ref. j65[. 

The outline of the remainder of this paper is as follows. In Sec. II we give the basic 
ingredients needed for our computations, including the weak Hamiltonian (Sec. II A) and 
expressions for the hadronic tensors and decay rates (Sec. II B). In Sec. II C we give a 
detailed discussion of the endpoint kinematics and light-cone variables, and in Sec. II D we 
briefly summarize a few results obtained using the optical theorem for the forward scattering 
amplitude, and the procedure for switching between partonic and hadronic variables that is 
relevant in the endpoint region. In Sec. Ill we turn to the discussion of the SCET formalism. 
Many of the ingredients necessary for our computation are readily available in the literature. 
Of particular note are expressions for the heavy-to-light currents at (9(A^) [66], which we 
have verified. In Sec. IV we review the derivation of the factorization theorem at LO, 
but do so in a way that makes the extension beyond LO more accessible. We consider 
power corrections of 0{X) in Sec. V, and show that they vanish. In Sec. VI we discuss the 
true NLO factorization theorem, which is (9(A^). In particular, in Sec. VI A we switch to 
hadronic variables and re-expand the LO result, in Sec. VI B we enumerate all the time- 
ordered products that occur at this order, and in Sec. VI C we show that the tree-level 
matching is simplified by using the SCET formalism. In Sec. VI D we give definitions 
for the non-perturbative shape functions that appear, and then in Sec. VI E we derive 
the factorization theorems for the most important contributions in some detail. Finally, 
in Sec. VI F we summarize the hard coefficient functions for the subleading time-ordered 
products. Next, in Sec. VII we present a useful summary of the NLO decay-rate results, 
including the phase-space corrections. We also compare with results in the literature where 
they are available. Our conclusions and discussion are given in Sec. VIII. Further details are 
relegated to appendices, including the expansion of the heavy quark field and calculation 
of the power-suppressed heavy-to-light currents at tree level in Appendix A, and a review 
of constraints on the currents from reparameterization invariance in Appendix B. For the 
reader interested in getting an overview of our results while skipping the details, we suggest 
reading Sees. II, III and IV, the introduction to Sec. VI, and Sees. VI A, VI B and VII. A 
reader interested only in final results may skip directly to the summary in Sec. VII. 
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II. BASIC INGREDIENTS 

In this section we give the ingredients necessary for studying the decays B —>■ Xg'-f and 
B Xud-v in the endpoint region to NLO. A proper treatment requires a separation of 
the scales ^ ml ^ m^AQCD ^ ^qcd the form of a factorization theorem. This is 
accomphshed by the following steps: 

1) Match on to the weak Hamiltonian, H\y, at /i^ = and run down to /x^ = m^, just 
as in the standard local OPE. 

2) Match Hy^ at jj? ^ on to SCET, with coUinear and usoft degrees of freedom and 
an expansion in A = a/Aqcd/'^6- Run from /x^ = to /i^ = /Xq ~ mbAqcD- 

3) At = /Iq integrate out the collinear modes, which, given a complete factorization in 
step 2), is trivial. Then run from /x^ = /ig to /i^ ~ 1 GeV^ ^ ^qcd- 

In Sec. Ill Al we discuss the weak Hamiltonian. The kinematics and differential decay rates for 
the endpoint region are given in Sec. Ill (^1 Then in Sec. II I II we give the necessary effective- 
theory Lagrangians and currents to 0{\^). 



A. Weak Effective Hamiltonians 



For B — > Xy^lu the effective Hamiltonian is simply 



eff 



V2 



(10) 



and the current w^^Pl}) is the basis for our analysis of the QCD part of the problem. The 
Hamiltonian for the weak radiative decay B —>■ Xgj is more involved: 



K 



4Gf 



eflf 



V2 



(11) 



1=1 



where Pr^l = 1(1 ± 75) and 



07, 

0{3,5} 



mbsaauF'^'' Pnb 



167r2-^"^"'^^^ 
{c^^PLh){dYPLc), 
id^,PLb)J2iqrPL,Rq) 



m,sa.,G^''PRb 



02 = {c^^.PLb){sYPLc), 
C?{4,6} = {df37''PMj2{qarPL,Rq(^) 



(12) 



with F^i, and G^^^ the electromagnetic and QCD field strengths. In the expression for TCl^, 
we have used unitarity, V^iy*^ = —V^f^V^l — V^^KTs^ have dropped the small (~ 2%) 
corrections from V^iy*^ and ms/mf,. 

For the total B — >■ Xg'y rate the perturbative corrections are known at NLO 

SQ- Ef- 
fective scheme-independent coefficients Gj% are defined in a way that includes contributions 
from the penguin operators (Ca-e). A totally inclusive analysis is considerably simplified by 
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the fact that at leading order in l/m{, the matrix elements can be evaluated directly in full 
QCD rather than first having to match on to HQET. For an endpoint analysis, the matching 
at uiw and running to /i ~ m;, is the same. However, at the scale /z ~ m;, the operators 
in Hw need to be matched on to operators in SCET before the OPE is performed. In per- 
forming the matching, the only subtle complication is the treatment of the charm mass. For 
simplicity, the approach we take here is formally to let rric ~ wij,, so that charm-mass effects 
are all hard and are integrated out in matching on to SCET. This agrees with the treatment 
of the Oj^y advocated in Ref. [g^ for the endpoint region. Since numerically ~ m^A, 
perhaps a better alternative would be to keep charm-mass effects in the operators of SCET 
until below the jet scale m^A. This second approach is more involved, and in particular it is 
clear from Ref. that it would necessitate introducing two types of collinear charm quark, 
as well as soft and ultrasoft charm quarks. For this reason, we stick to the former approach 
and leave the latter for future investigation. 

At lowest order in the A power expansion, there is only the SCET analog of the sa^i^Pub 
current called J(°) (cf. Eq. dMD), and O i_8 can make contributions to its Wilson coefficient. 
At NLL order in the effect of the other operators can be taken into account by using (6^ 

3 



c,i,, z) = cr\,) + i^cr\,) + E cr\,) [nip) + in (^)] , (13) 

in place of C7 when matching on to SCET. Here the dependence on p = ml/ml enters from 
the four-quark operators with charm quarks. In the endpoint region the 0{as) effects from 
the process b s'jg all appear in the jet and shape functions. For later convenience, we 
define 



C7 [p) 



The photon in B ^ Xg'-f is collinear in the opposite direction to the jet Xg, so propaga- 
tors connecting the two are hard. Thus, beyond LO in the power expansion, the photon will 
typically be emitted by effective-theory currents J^'^^ (which could be four-quark operators). 
We shall discuss the matching on to these subleading currents for O7 only. Some of the 
contributions from the other Oi will just change the Wilson coefficients of the subleading 
currents and thus not modify the structure of the power-suppressed factorization theorems 
(indeed some of them are already known since they are fixed by reparameterization invari- 
ance). These other operators may also induce time-ordered products that would involve 
operators with quarks collinear to the photon direction, but these are not considered here. 



B. Hadronic Tensors and Decay Rates 

In this subsection, we summarize general results for the hadronic tensors and decay 
rates, without restricting ourselves to the endpoint region. For both decays B Xgj and 
B — ^ Xuii^, momentum conservation for the hadrons gives 

p^ = mB^;^=p^ + g^ (15) 
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where Px is the sum of the four- momenta of all the hadrons in X, is the momentum 
of the 7 or the pair of leptons (^z^), and the velocity satisfies = 1. For B Xs'j, 
gfj. — 7^. g 7^/^/2 = E^n^, where v ■ q = is the photon energy in the rest frame of the B 
meson. In this case = and 



mB{2Ex, - ms) = rriBimB - 2E^) , (16) 
so the differential rate involves only one variable, mx or E^. For B X„£z/, Eq. (jl5|) implies 

(17) 



Ex = v-Px 



2m 



V q = niB - Ex = Ei + E^ 



B 



where pj^ = m?x. The differential decay rate involves three variables, several common choices 
of which are {E^, E^, q"^}, {Ei,v-q, g^}, or {Ei, m\, g^}. 

To derive the inclusive decay rates for B —>■ Xgj and B —>■ Xuii^, the matrix elements 
are separated into a leptonic/photonic part and a hadronic part VF^y. Here 



= ^J2(^^^'^\PB-q-Px){B\Jl\X){X\J,\B) 



X 



g/MuWi + Vf,VuW2 + ie^uaf^v^q'^Ws + q^^quWi + {v^q^ + Vyq^)W^ 



;i8) 



in which we use the hadronic current J and relativistic normalization for the \B) states. For 
convenience we define projection tensors P^^ so that 



They are 



1 /lu v^v^ + q^q^ — v-q{v^q'' + v^q^) 
'2^ ^ 2[g2- (t;.g)2] 



(19) 



(20) 



-'2 



-^4 



3g2 ^ + q^g^"' - q'^q'' 
[g2 - [v-qf] 

[g2-(i;.g)2] 



pfiu 
^3 



-ie'""''^qaVf3 
2[g2 - [vqy] 



q'pr 



2v-q 



Contracting the lepton/photon tensor L^^ with W^^'^ and neglecting the mass of the 
leptons gives the differential decay rates 



8E^ 



dE^ 



ru 
n 



nil 



96 



(4iyf - - 2E^Wr 



(21) 



B 



q'W^ + {2E,E, - qy2)W^ + q\E, - E,)W^ d{AE,E, - q^ 



dE^dq^dEy 

where Wi = Wi{q'^,v-q) and the normalization factors are 



327r4 



VtbV^ll'^ ttem [mb{mb)]'^\Cf''^^ {mb) \ 



^ 



1927r3 



ub 



(22) 
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Here mb{jj) is the MS mass. For convenience, we have pulled out a Wilson coefficient C^^^'^^ 
so that contributions from other coefficients appear in ratios Q/Cf in the SCET Wilson 
coefficients (for example the quantity in Eq. ()14|) ). Note that we have chosen to stick 
with hadronic variables here (using uib rather than mf,). When we eventually compute the 
Wj, we shall have to deal with switching between partonic and hadronic variables. However, 
we shall see that the situation is quite different from that in the local OPE (as we discuss 
further in Sec. Ill Dl below). In particular, it is the hadronic phase space that turns out to 
be required. 

In Eq. (1^ . < El, Ey < {rn\ — m^)/(2mB). A set of useful dimensionless hadronic 
variables is 



7 2^7 g2 ^2 

Xh = , xh = , Vh = , sh = ■ (23) 

niB ttlb mj 



In terms of these variables. 



2E,j 2Ex , \ 

I - Sh + yH - xh , = 1 + Sh - Vh , (24) 



niB ruB 
and Wi = Wiiyn, sh)- For B X,7 

dT' 2a;, 



^{4W^i - - x\Wl)^ , (25) 



with 0<x]^<l- m\,/m%. For B XJu, 

, "^'^^ = T^o'^AmB\yHWl' + l[{l-XH){xH-yH)-XHSH]W^ (26) 
axH dyn dsn 12'- 

+ ^mByH{2XH + SH-yH-'^)W^'^9[{l-XH){xH-yH)-XHSH] , 

and, depending on the order of integration, there are several useful combinations of the 
limits, which are shown in Table HI If we integrate over all values of xh (cases iv) & v)), the 
rate becomes 

f^^ = 2mB^{l-yH + SH?-^SH \l2yHW^ (27) 
dyH dsH y 

+ [{l-yH + SHf-ASH]W^]. 



C. Light-cone Hadronic Variables and Endpoint Kinematics 



We are interested in the jet-like region corresponding to Aqcd ^ f^x <C Ex for both 
B and B — >■ X^iu. In this region, the hadrons in the X occur in a jet in the B rest 

frame with Ex ^ mB and mj^ ^ tubAqcd-^ The momentum of the states X is therefore 



^ In Refs. |45l l46l| an intermediate situation, < with ^msAgcD <C A <C ms, was also considered, 
but we do not consider it here. 



11 



i) 0<XH<l-rl 

a) < XH < I - rl 

in) r"^ < sh < ^ 

iv) r1 < sh 

v) Q<yH<{l-r^f 

vi) 0<yH < (i-'^tt)^ 

TABLE I: Limits for different orders of integration in B ^ Xu^u with variables {xjy, sjy, yjy}. 
Here = m^/niB, while = [(1 + yn - sh) ± + yn - snf - 4yH]/2 and 

l^max* ,j^mm*| _ j^-max^ ^mmji^^^^^ . Results for the phase-space limits of partonic variables are 
obtained by dropping the i7-subscripts and setting r^^ = 0. 

restricted, but they still form a complete set for Eq. ()18p. The width of the jet is determined 
by noting that the typical perpendicular momentum between any two final-state hadrons is 
Ap±_ < y/niB^ ~ 1.6 GeV, where we use A ~ 0.5 GeV to denote a typical hadronic scale for 
i?-mesons (examples being Aqcd and A). We also assume that there are enough states with 
~ m^A that the endpoint region is still inclusive and Eq. is valid. As we shall see 
below, the factorization in this region gives non-perturbative shape functions rather than 
just local operators. 

For jet-like kinematics it is natural to introduce a light-like vector n along the jet direction 
and an auxiliary vector n such that = fi^ = and n ■ n = 2, and refer to components 
= {n-p,n-p,p'^). Some useful decompositions are 

9T = 9"' - lin^^n^ + 9^" = 9'" ~ v^v^, 

pi^ = p^' - !1J^ n^" - ^ , p>^ = p^" - v^v-p , 

where we take e°^^^ = 1. Note that the subscript T means transverse to v'^, so pt^ P±- For 
the final factorization theorem for the differential decay rates we shall use a frame where 
q^ = v^ = 0, and = {n^ + n^^)/2.^ Thus = n-qn^ /2 + n-qn^^ /2 and 

n-q = niB — n-px , n-q = uib — n-px ■ (29) 

For B Xgj the photon momentum is taken along the n light-like direction, i.e. q^ = E^n^, 
and 

n-px = mB, n ■ px = mB ~ = mB{l - xIj) . (30) 



^ If one desires, he or she can take v = (1,0,0,0), n = (1,0,0,-1), and n — (1,0,0,1). A more general 
frame is required only for working out the constraints from reparametcrization invariance. 



rl<SH<l-XH ^<yH<XH- 

^<yH<XH-i^ rl<SH<^{xH-yH) 

0<XH<l-SH 0<yH<XH- 1^ 

0<yH<il-V^f xf''<XH<xf^- 



rl<SH<{l-^f xf'' <XH <x'^-- 



H 

H <XH<x'ft-* rl<SH<l + yH-ff--XH 



(28) 
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i) 
''J 


< xu < I — r"^ 


T"^ liu 1 — T,u 
' 71" — — ^ •^H 


mcLX-|^ 1 — 




UH > — 


< 1 


ii) 


< XFf < 1 - 


inaX'fl — xu T—-^ — \ ^ U TT ^ \ 


— < 

Vh - ^ 


< 1 


— Xu 




iii) 




< xh < 1 — uh 


maxjl — 




vi;] < Vh 


< 1 


iv) 


rl<UH <l 


2 






<1-uh 




V) 






i-yH< 


XH 


<l-UH 




vi) 


<yH 


l-yH<^H<l-^ 




< 1 


- XH 





TABLE II: Full phase-space limits for B X^H-v with variables xh^ Vh^ uh- The parameter 
r-jr = m-ir/n^B- Results for the phase-space limits of partonic variables are obtained by dropping 
the i7-subscripts and setting = 0. 

For B Xuiu the phase space is more complicated and for convenience we define the 
dimensionless variables 

Vh = 5 Uh = . (31) 

rriB rriB 

Now = n-pxn-px and n-px + n-px = {f^B ~ <f + fT^\)/^B-, so 

SH = UHyH, Vh = -uh){1 -Vh) , (32) 
and, making the choice IJh > uh, we have 

[yn, Uh} = ^ 1 - yn + sh ± a/(1 - y// + sh^ - 4:Sh ■ (33) 

So far we have not made any restriction to the endpoint. The variables y^j and uh provide 
an equally good description of the full B X^d-v phase space as yn and sh, namely 

1 d^F" f 1 

+ ^{l-u„)(l-y„)(2x„ + u„ + y„-2)W;Y (34) 

where Wi = Wi^UHyVn) "we have suppressed the theta function from Eq. (j^ . Integrat- 
ing over Xh from Table ITTl (cases iv) & v)) gives 

fi'W^ = 2AmB{yH-UH?{{l-UH){l-yH)W^ + ^{yH-uufW^] . (35) 

The full limits for y^ and uh are given in Table IHl 

In Ref. jHl, it was pointed out that a natural set of variables in the endpoint region 
consists of the hadronic variable n ■ px and partonic variable n ■ p, where 

n-p = n-px + nib — uib = n-px — A -|- . . . , 



n 



P^ = -jn-p+ —n-px- (36) 
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m. 



SCET 
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1 r. 



.6 ■ 



.4 ■ 



annihilation 



U 



H 



charm .6 
contamination 



charm 
contamination 



annihilation 



SCET 
region 



0.2 0.4 0.6 0. 



THj^ (GeV) 



FIG. 2: Allowed phase space for B — > X^^v, where < mx < rns- The second figure shows 
the same regions using the variables defined in Eq. pij) . We indicate the region where charm 
contamination enters, mx > mu, and the region of phase space where annihilation contributions 
enter. Also shown is the region where the SCET expansion converges, which is taken to be uh/vh — 
0.2 and corresponds to m\/{AE'jr) < 0.14. 

They are natural because the LO factorization theorem dictates that the kinematic variables 
appear in the jet functions and soft functions only as 

J{n-pk+), fXn-px-k+), (37) 

where is the convolution parameter (cf. Sec. II V|) . We shall see that this remains true 
of the shape functions and jet functions at subleading order in the power expansion. For a 
dimensionless version of n-p we use y = n-p/nib- A comparison of the phase space with the 
variables {yn = Q'^/m'^, ^/sh = mx/rriB} and {uH,yH} is shown in Fig. Eland corresponds 
to the limits shown in columns iv) and v) of Tables H] and El The figure on the right is the 
analog of Fig. 1 in Ref. (7l| with dimensionless variables. 
For a strict SCET expansion we want 

n-px 



< < 1 



(38) 

n-px 

where the expansion is in the parameter X^- Eq. (jHH)) is equivalent to Uh/Vh ^ -^h- Fo'^ 
B Xs'j, Eqs. ()16|) and ()38|) imply that the endpoint region is 



For B X, 



> (mB/2 - A) 
satisfying the criterion in Eq. 

mx < Xh 



2.1 GeV. 



with values ~ 



2E 



X 



1 + A2 



0.37. 



(39) 

0.2 is equivalent to 

(40) 



H 



OT yn ^ I — 2.Q8y/sH + sh- Not exceeding a given expansion parameter in Eq. corre- 
sponds to specifying a triangular region of phase space (shown for uhIVh — 0-2 in Fig. 
We refer to this as the SCET region of phase space.'' As can be seen from Fig.|2l the simpler 



^ Note that the SCET expansion here is actually in powers of since odd terms in SCETi tend to be 
absent [t^. In Sec. |3 we show that the C(A) contributions for inclusive decays indeed vanish. Thus 
— 0.2 is not a large expansion parameter. A smaller value Xjj — 0.1 could be chosen if desired. 
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restriction yn ^ ^ — 2.44^^5/7 gives a very good approximation to the SCET region since 
the boundary is roughly hnear in the and mx variables. 

In calculating decay rates at subleading order, it is important to define carefully how the 
phase-space integrals are treated once we compute doubly differential or singly differential 
decay rates. The philosophy we adopt for B Xudv is that we use SCET to compute the 
Wi, and hence the triply differential decay rate, for the SCET region in Eq. ()38|) . In general, 
one may wish to integrate this rate over a larger region of phase space, and thus need to 
construct the full VFj's. This could be done using 

^fuii ^ wr^^e{X'^y^ - un) + Wr^'eiuH - Xlvn) , (41) 

where the SCET expansion in A is used for the first term and the standard local OPE in 
A/mfe is used for the second term. Thus, \h does not play the role of a strict expansion 
parameter, but rather provides us with a means of interpolating any differential spectrum 
between the full OPE and full SCET results by varying \h between and 1. We only 
consider the first term in Eq. (PT|) here. Depending on the final spectrum that one looks 
at and the other cuts imposed, the error in including a larger region of phase space than 
the SCET region may be power suppressed. The parameter provides us with a way of 
testing this by considering the difference between taking Xh = 0.2 and \h = 1- We present 
our final results in a manner that makes it easy to take the Xh — ^ 1 limit for situations 
where a large enough region has been smeared over that this is the case. 

One can also refer to a shape-function region, corresponding to the region < uh ^0.1 
where the non-perturbative function / is important. The expansion for B Xg'-y to second 
order in A, where 1 — xjj ~ A^, is 



dx]j ° niB 



- - TUB } - {l-xl)^AW^-W^-2mB W^^ 



(42) 



For B Xuity, if we integrate over all xh and expand in uh ~ A^, then from Eq. the 
first two orders in the expansion are 



24:171 B 



dyjj duH 



4 ^ ^ ttS 



(43) 



By the endpoint region in Xh we mean x'jj < Xh < 1 — r^, where 1 — x^j ~ 0.1 corresponds 
to making a cut on the lepton's energy spectrum. The limit on uh forces it to be small, 
Uh < ^~^'hi so shape-function effects are important here. This cut still allows a large range 
for yjj. Expanding the triply differential rate in Eq. ()34j) and keeping the first two orders in 
the expansion for 1 — xh ^ uh ~ A^ gives 



Fq dx^ dyn dun 



24m 



B 



{yH{i-yH){w^ + '^w^)} (44) 



-UH[{i-yl)w^ + '^yni'^-yH-yDw^] 
2 

+ [1-xh-uh){^-W^ -mByH{l-yH)W^] 
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Finally, to consider the dV /{dq^dw^x) spectrum in the endpoint region we let = l—yH + sh 
and expand in sh/C'^, which to linear order gives 



{i2C(i-c)w^r + cX"} + ^{2(c'+2c-2)w^r - c'w^} ■ 

(45) 

In each of Eqs. ()42|1 - ()45|) there will also be an expansion of the Wi themselves, which we 
discuss later on. 

The results in this section can easily be extended to any desired order in A /tub by 
expanding to higher order in the phase space. 



dyn dsH 



2m 



B 



D. OPE and Partonic Variables 

We have not yet made use of quark-hadron duality or formulated the method for comput- 
ing the Wi. The usual procedure to compute the Wi is to use an operator product expansion 
and calculate the forward scattering amplitude 

where 



T 



-I y"rf^xe-'«-"Tjt(x)J,(0), (47) 

and for the Tj^^'^\q^,v-q) we use the corresponding hadronic currents, which are 

j; = sta^^q'^Pn b, j; = u YPLb . (48) 

Here comes from the operator O-j. The operator product expansion relates the W^^'^^ to 
the forward scattering amplitudes through 

Wi = --lraTi. (49) 

TT 

When we compute the Wi with an OPE, the partonic variables depending on and the 
hadronic variables involving will need to be related order by order in the l/m;, expansion. 
In particular, the heavy meson mass to second order is 

„, = ,„. + A-_^-!££:MMl) + ..., (50) 

where we shall also use A2 = Ci;'(/i) A2(yu) as a definition of the non-perturbative matrix 
element that has the perturbative coefficient cpifi) absorbed. 

In applying the local OPE, part of the expansion involved in switching to hadronic vari- 
ables occurs because the phase-space limits are partonic. In fact, if we calculate the triply 
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differential rate for B Xuijy with the local OPE, and then consider integrating it over the 
hadronic phase space, then the integrand has support over the partonic phase space only, 
so the limits are reduced to this more restricted case. In the local OPE the signal for this is 
the occurrence of factors like 

S^-\{m,v - qf] , (51) 

which must be smeared sufficiently so that quark-hadron duality can be used. For the 
6[{mbV — g)^] that occurs at LO, integrating once to get the doubly differential rate gives a 
theta function that imposes partonic limits. 

On the other hand, with the SCET expansion in the endpoint region the support of 
the triply differential rate is larger. We never encounter singular distributions like the 
one in Eq. ()51|) . but instead obtain a non-trivial forward 5-hadronic matrix element that 
gives f^'^\i~^)- This function knows about the difference between the hadronic and partonic 
phase space already at leading order in the power counting, and more generally the LO 
factorization result with 0{as) corrections (cf. Eq. below) does not cause a restriction 
of the hadronic phase space. Therefore, we shall use the full hadronic phase-space limits in 
our computation.^ 



III. SCET INGREDIENTS 

To derive the expansion of the Tj in the endpoint region, we determine the SCET currents 
and Lagrangians in a power expansion in A, and use this expansion to separate the collinear 
jet-like effects from the non-perturbative ultrasoft shape functions. This will allow us to 
determine results for the Tj order by order in the expansion. We write 

T, = 7;^°^ + t/') + Tf ) + . . . . (52) 

The ingredients required for this include expansions of the Lagrangians and currents to 
second order in A, 

c = + + £(2) , J = + + , (53) 

and are described in this section. 

The LO Lagrangian for usoft light quarks and usoft gluons, CuJ , is identical to full 
QCD. For usoft heavy quarks we have the HQET Lagrangian, which has an expansion in 
A /nib ~ A^, namely 

C'^^^ = hyiv-DusK, (54) 

1 1 
C'f^^ = Oh = hv{iDTfhj, + -CF{fi)hv<y,iiy gG'^"^ , 

® Note that this might imply that direct calculations of less differential subleading spectra must be treated 
with care. For example, if one directly computes a singly differential rate by tying up lepton lines then 
the partonic phase-space restrictions might appear to creep back in if one is not sufficiently careful about 
the structure of the factorization theorem. 
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where Ci?(/i) is the Wilson coefficient of the chromomagnetic operator. At leading order, the 
SCETi Lagrangian for coUinear quarks is ji^ li^l 



(0) 

ii 



in 



in-D + ipj:W^WUpt 



in ■ 



(55) 



where the collinear covariant derivative is iDf! = + gA'^, with label operator V^, the full 
derivative is in-D = in-d + gn-Aus + gn-An, and the Wilson line W is built out of n- 
fields: 



W 



perms 



5Z {- ■^^■^n,g{x) 



(56) 



The collinear covariant derivative satisfies f{in-Dc) = Wf{V)W'^. In the leading-order 
collinear Lagrangians for both quarks and gluons, C^^^ + C^g , we can factorize the usoft 
gluons with the field redefinitions |37| 



A^„ YAtY^ 



(57) 



Y = Y [x) = P eicg iig I ds n- Aus{x+ns] 

J —oo 



where the new collinear fields no longer transform under usoft gauge transformations. For 
convenience, we also define 



H = Y^'h 

The leading-order collinear Lagrangians become 



Vr = W^DrW . 



Y^D.,,Y 

^lw^[tn-D,,tD^]wY (58) 



C 



(0) 



Xn 



2 



4? = -;jtr[GrG^g+g.f., (59) 



where igG'^'^ = [iD'^^iD^] and g.f. stands for gauge-fixing terms. Thus the field redefinition 
removes the coupling of n-Aus gluons to collinear fields at LO. At subleading order it 
organizes the usoft gluons into gauge invariant blocks, like (y^...y), which are inserted 
only a finite number of times. The subleading Lagrangians that we require are all available 
in the literature |73j,i7^]. After one makes the above field redefinition, the subleading quark 
Lagrangians we require are j?^ 



Hi 

Hi 



Xntrts^^Vc^Xn + h.C. 

Xn igl$^i)us + h.c. , 
1 ^ 



(60) 



XnlKsWus p^Xn, 

xJ'Ptin-Vus:^^ | ift Xr. 
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We also need the subleading terms in the mixed usoft-coUinear gluon action, 

= -2tr[^r^;]^?j-, + g.f., (61) 
^(2a) ^ [u^'^u^^] gf^^ gl - tr {W^n;-] gt - tr [^r^r] 9^ 9^r + S-f • , 
4f = -tr[^^7^;]r^,n. + g.f., 

where 

ig = [wUD^W, YHD-^^Y] , (62) 

and g.f. denotes gauge- fixing terms that are required by reparameterization invariance. 

Next consider the heavy-to-hght currents. For simphcity, we work in a frame where 
v± = 0. Expanding the J" current to (9(A^) gives 

+ J2ldcoAf\u)jf\co)Y (63) 

where the superscript (0), (1), (2) indicates the order in A. We have an analogous result 
for with (v) — >■ (t) and m — ^ s. After one makes the field redefinition in Eq. (|57|). the 
leading-order SCET heavy-to-light current is js^ 



Jf\uj) = {^nW)^rf\Y^K) = Xn,.Tf'H. , (64) 

where j = 1-3 for u and j = 1-4 for s. Despite the fact that the operator in J^^^ is 0{\^), we 
use the superscript (0), to indicate that it is the lowest-order current. At the first subleading 
order, we have the C(A^) currents 7^ 74. TgI 77| 

^r(^) = l^{xniD:^)rf'''>n., (65) 

J-^^\uJi, UJ2) = — ^ — Xn,c.i {igB"j_)^^ rf^^ Hy . 

The subscript u notation in Eqs. (jMj) and indicates that the field carries momentum 
u;, for example 

Xn,u = Xn 5(^ - P^) , {igB'^cA.U = iwB^J 5(^2 " P^) • (66) 
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Finally, at second order we find that the basis of 0{X^) currents consists of 



2mb 
1 



Xn,uj ^ja ''^ us '^ti 



n-v ign-Bc 

Xn' 



V n-v 
-1 - ( ign-Bc 
mfe V n-v / u)2 
-1 



.{2c) 
j 



nibn-vui 
-1 / 1 



= -Xn,.rf^^ ^Vi-n., (67) 

j{2d) , 

J] '[yji,uj2 

These currents were first derived in Ref. j66| (see also [73|) and we agree with their results. 
Two differences are that in Eq. (j67j) we have determined the functional dependence and 
the number of parameters tUj, and that terms from the separation of label and residual 
momenta (or equivalent ly the multipole expansion of the fields) are treated differently here 
from Ref. 6^| . In Appendix El we give a more detailed comparison, along with the details 
of our calculation. 

For the LO Dirac structures we use the basis from Refs. 74, 7^, which is the most 
convenient for considering reparameterization invariance, i.e. 



Xn,ijJi 1 ~ - 

mb \n-vV 



p«(0) 
^ 1-3 



n"^ 



n-v 



n-v n-v 



(68) 



Here we slightly reorganize the basis in Ref. [t^ to reflect the constraints from RPI and to 
use subscripts j > 11 for the currents that vanish when v± = 0. This gives 



l^f^A ^,ltJ.,„M^,a ^,OilfJ-,^,iy] ^°^[M.„,!^1 



(69) 



Note that we do not list evanescent Dirac structures that can become necessary when com- 
puting perturbative corrections in dimensional regularization (see for example Refs. 7^ 79|). 
For our purposes, the Dirac structures for the J^^-* currents that appear at lowest order in 
as(mfe) will suffice; these are 

i 1 —^B.l —^Rl±lpi —^Rl, JL 1 — -r/J 7 — , 



rf'^=PRr7H 



^s{2a) _ p ■UT.n y.s(2fe) _ p „ ^ „-^A«r^ nv,s(2c) _ p • ur^ y.s(2(i) _ p ■ ur^ 4^ 

i 1 — -TR ta"^ 'yrpqT , 1 1 —'t^Rl±2 ' ^1 —^R It 1 J- 1 — -tr icr"^ 



^s(2e) _ p UT i \ ^^(2/)_p • ur I I 

-■-1 — ^Rla 2 2 7/3 5 J- 1 —J^^RT'O Qr'Ja 7/3 



(70) 
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The complete tree-level currents for arbitrary v± are given in Appendix^ and in Appendix iBl 
we show how the v± dependence is necessary to satisfy the full constraints from RPI at this 
order. When the 0{\^) currents that show up at order as{m}j) are determined, the number 

('2') /fTTTTl* 

of possible structures will increase, as occurred in Eq. (jH^I)- However, this will not 
affect the form of the factorization theorems we derive, since they are expressed in terms of 
traces of the T^- factors; this makes it trivial to incorporate new Dirac structures that arise 
beyond tree level. 

At lowest order in a<((m;,), the non-zero coefficients in Eq. are 

Cf^) = C? = 5f ^ = 5g = -eJ") = 1 . (71) 



The one- loop results and RGE-improved coefficients can be found in Ref. [39|. Our 
coefficients are linear combinations of these: 

C!")(^.l) = l_^£(!^|21„V) + 2Li,(l-^) + l„M(ff^)+g + 6}, 

O2 [UJ, ij 



An [{l-uj) 
as(mfe)C'ir f (1 — 2cj)u) ln(u)) lj 



78| . and the anomalous 

m. 



Cf\u;,l) = l + A,(m,,p)-^i^^^(21n2(c^)+2Li2(l-^)+ln(c2;)(^ . . , 
C?(^,1) = 0, C1V,1) = ^^{^^}. Cl"(^,l) = 0. (72) 

where lj = 00/ nib and A^(mfe,p) was given in Eq. Thus, non-zero values for €2^^ and 

Cf^ are generated at one-loop order, while 6*2*4 are still zero at this order. For the J'-^*^ 
currents, the one-loop matching coefficients were derived in Ref. 
dimensions and RGE-improved coefficients were computed in Ref. 

By reparameterization invariance b[^\{uj) = C[^\{uj) and Bi}_^{uj) = C[^}^^{uj), so their 
one-loop matching coefficients are in Eq. (f7^ . For the (9(A^) current coefficients the RPI 
constraints are 

= A? = -^4^ = CfV) , 4^ = 4^ = 4^ = 5?Vi,^2) . (73) 

Note that we have not included effects associated with a non-zero strange-quark mass in 
our basis of subleading SCET operators, and it would be interesting from a formal view-point 
to consider the possibility of such power-suppressed terms in the future (cf. fiH). 



IV. LEADING-ORDER FACTORIZATION 

In this section we review the leading-order factorization theorem [s^ as derived using 
SCET js^. The result is discussed in detail, which will allow us to refer back to this section 
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„ , „ , „. Hard, Jet, and _ ^ 

i -product Example Diagram Usoit Operator 

Shape Functions 



^ 

TABLE III: Lowest-order insertion of SCET currents. The double hues are heavy quarks and the 
dashed hne is a cohinear hght quark. 

when some of the steps are repeated at NLO. Throughout this section we shall suppress the 
u and s superscripts, since all the manipulations are identical in both cases, B — > X^ii^ and 

To derive the leading order t/^^ we need only consider T)i°\ given by taking two LO 
currents J^^^ , as shown in Table 11111 The phase factor in Eq. (|63|) gives 

and exp(— iP^ ■ x). The current has q± = 0, so this V± term will contribute only to fixing 
the perpendicular momentum of the jet function to be zero. The large label momentum in 
the phase in Eq. ()63|) also gets fixed by momentum conservation: 

V = rrib — n ■ q = nib — f^B + n ■ px = n ■ p , (75) 

where n-px is the large momentum in the jet X. Then the remaining momentum r'^ ~ A 
since n ■ r = n-q + V ~ nib = and 

n-r = n-q ~ iTLb = niB — nib — n-px ■ (76) 

At lowest order 

n-p = n-px, n-r = A — n-px, (77) 

where both A, n ■ px ~ A (and higher-order terms in — nib will be needed only when we 
go beyond LO). For the time being we stick to the partonic variables n-p and n-r; later, we 
shall perform the expansion involved in switching to hadronic variables. Using the states 
defined with HQET, we get 

WlH^ = ^Imi(5.|fW|5.), (78) 

f(°) = -^jd'xe—TJf^{x)Jf}{^). 

Separating out the hard Wilson coefficients, we have 

fl^J = Y. jdujduj' Cy{uj')C,{uj)5{uj'-n-p){-t) jd^xe-^'-^ T jf^{uj\x) jj2(a;,0). (79) 
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The effective-theory currents in the remaining time-ordered product depend only on colhnear 
and usoft fields describing momenta <^ m^, i.e. 

T jf^ [J, X) jf} (a;, 0) = [njf^Xn,^.'] (X) [Xn,a;rfjK] (0) , (80) 

where f = 7°r^7'^. It is useful to group the collinear and usoft fields into common brackets 
by using a Fierz rearrangement. For spin and color we can use 

1 ^ 

l®l = -Y,F^®Fk (81) 

fe=i 

2 L4iV/ 4^ ^ 2NJ^ 2 ^ ^ 2N, 2 ^ 

+ (^T")®(^) + (-^7,r^)®(^) + (-^7!T'^)®(^) 

Equation (j81|) is valid as long as the identity matrices on the LHS are inserted such that on 
the RHS the FJ^ appear as ^n^FJ^BC,n, where A and B do not contain ^ factors. In Eq. (j8(jp 
we insert identity matrices to the right of the F^?^ and to the left of the ^fj, which gives 

{B^\TjfJ{u;',x)J^'J{u;,0)\B^) = (5.|T[7i.(x)rJ?>," Ffj 7^,(0)] 

x{0\T[xnAO)F^Xn,..'{x)]\0). (82) 

Here the states \B^) have HQET normalization, {B^{k')\B^{k)) = 2v^{27rf5^{k'-k), and 
are defined as energy eigenstates of the LO usoft Hamiltonian generated from Cus- 
Only the vacuum matrix element with A; = 1 is non-zero:^ 

{0\T[xnA0)^Xn,.'ix)]\0) = {0\T[i^^WU0)^^iW^^n)^ix)]\0) 
= (-2.) 6{uj~uj')6'{x^)6{x-') l^e-^'^^-/' ji'\k^) . (83) 



This definition of J'^^'^ agrees with Ref. j33|.^° Owing to the forward matrix element, the 
momentum-conserving delta function gives the S{uj' — u), which can be used to eliminate 
the Lj' integral in Eq. (fTSj). and the delta function in Eq. (f7^ then sets 

u = n-p . (84) 

The appearance of the 6{x^) ensures that in the remaining usoft matrix element in Eq. ()82|) 
time ordering is the same as path ordering along x~ . Up to order a^, the one-loop diagrams 
plus counterterms give 



a^(/i)CV 



,(85) 



9 Note that in Eq. JHSJ 6^{xj_) ^ A^, since from the momentum-space view-point it is obtained by combining 
integrals over both the label momentum p± ^ X and residual momenta k± ^ X^, rather than just the 
latter. 

To check this, one must note that here we contract color indices on the LHS and have an extra minus 
sign, since the x„ fields are in the opposite order to Ref. ■ 
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*' X 
xV ■> -Qi 




FIG. 3: Momentum routing for leading-order insertion of currents. 



which agrees with Refs. |45|, |4f 

Next, we simphfy the spin structure of the B matrix element in Eq. 



using the formula 



PTP 



P.Tr 



1 



1 



m=l 



(86) 

= 

for P^ = (1 + i))/2 and = Py'jJ^'-fr^P^. For our LO matrix element the term 



which is valid between heavy quark fields. Here, Pf = P^, Pf = Pv/2, P2 = s 



vanishes and, taking into account the delta functions in Eq. 

= n-x /2 , 

namely 



, we have only a function of 

(87) 



-(0) i p(o) , 

2 



(p,|r[H,(x) f;; f T^^j nM\ \b.) = Tr 



2 ^'^ 2 . 



{By\K{S:)Y{S:,0)K{0)\B,) , (88) 



where Y(x,y) = Y(x)Y'^(y) = P exp yig ds n- Aus{sn/2)j . Combining the phases from 

Eqs. ()78|) and (jHHj) and noting that d'^x = {dx'^dx^d'^x±)/2, we see that the matrix element 
in Eq. ()82|) is now 



i Tr 



2 ^'^ 2 . 



dk' 



'dx 
Air 



-ix-{r++k+)/2 't{0)(i.+ 



J^"^ik+){B,\K{x)Y{x,0)KiO)\By). (89) 



After we pull out the large phases, the residual momenta are as shown in Fig. El where 
i+ = r+ + /c+. Equation (jH^ involves just the leading-order shape function 



dx 
An 



-ix-£+/2 



{B,\h,{S:)Y{i,0)K{0)\B,) 



{By\K5{t -in-D)h.,\By) . 



(90) 



Note that we have taken mf, ^ A, and that Eq. (j9(J|) depends only on the residual momentum 
and masses that are of order A. Since we are free to integrate by parts in the forward matrix 
element, it is evident that f^'^\i^) is real. Momentum conservation implies that the residual 
momentum of the 6-quark should not be larger than the residual mass A of the P-meson, so 
f{i'^) is non-zero only for < ttib — rrib = A. Inserting r"*" from Eq. (f7H|) . this implies that 
k~^ ^ P~x- the renormalized shape function depends only on ratios of the dimensional 

regularization scale n and momentum scales ~ A. 
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To derive the leading-order jet function we need the imaginary part of {—i) times the 
result in Eq. (jHIIjl . which is the imaginary part of J^°\k+) /(°)(^+). Although the B matrix 
element integrated over x~ is real, it does provide a theta function that restricts < 
and therefore cannot be ignored in taking the imaginary part. The function {k'^) has 
an imaginary part for fc"*" > only, so we have support for the imaginary part over a finite 
interval only. The crucial point is that these theta functions are generated by dynamics 
within the effective theory and occur even though we have taken the m?, — oo limit. For 
the leading-order jet function for B ^ we therefore have [uj = n-p] 

^ lm\j^^\k+) e{p\ - k+)d{k^ 



J^'\uk^) 



(91) 



By RPI-III invariance m, this is a non-trivial function of only the invariant product 
{ujk~^), times an overall factor of u and the theta functions. For notational convenience, we 
suppress these extra dependences when writing the arguments of J'^^^ as w/c"'". 

Since r"*" = n-g — m?, < < A, or equivalently < A;"*" < pj, it is convenient to use the 
variables = n-p n-p x and where 



(92) 



so that ujk'^ = zn-pn-px 
taking the imaginary part of Eq. (jH3j) gives 



and J^u\k+) is non-zero only for < z < 1. At order a^. 



UJ 



pi { 

4:71 



1 + 



An 
A\nz 

z 



21n'4-31n4 + 7-7r2 



(93) 



4 In 



PZ 



1 



'i-z)e{z) 



with the standard definitions for the plus functions. Our definition for the jet function agrees 
with Refs. |45[l46| once we compensate for the extra 6'(pJ — A;+) that we included in Eq. ()9ip . 
The jet function depends only on the dimensionless parameter z and ratios of /i^ and the 
momentum at the jet scale, p^ ~ QA. 

At this point we have separated momentum scales in SCETj , and, combining Eqs. ()78p. 
fj82|l . (j83|l . (j9T|l . (j88|) . and (j90|) . we arrive at the LO factorization theorem in terms of partonic 
variables, 



(0) 



h,{n-p,m,,fi) / di+ J^'\n-p{i+-r+),fi) f<'\i+,fi) 



hi{n-p, mb,fi) 



dk+ J'^°\n-pk^,fi) f''\k++r+,fi) 



(0)/ 



Using Eq. ()77|) . the final LO result in terms of hadronic variables can be written as 



W 



(0) 



hi{p^,mi„n) 
hi{p^,mb,fi) 



4 



dk^ {p-^ k^,fi)r>{k^+A~p^, ^) 



(0)/ 







hi{px,mb,fi) p\ / dz J^^\zmj.,fi) f^\k - p\{l-z), jj) . 



(94) 



(95) 
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For practical applications we would also make use of a short- distance mass definition for mb 
(the cancellation of infrared renormalon ambiguities for the shape function was demonstrated 
recently in Ref. 0]). Here, the dependence on whether it is Xg or Xu occurs only in hf'^ 
and the values of the kinematic variables n ■ p and n ■ px- The functional forms of j'^^^ and 
/''°'* are independent of which process we consider. The hard coefficients are given by 

hi{uj,mb,n) = ^Cf{uJ,mb,^I)Cj{uJ,mb,^I) Tt Y^fl^'^f'^ ^^'^ ' (^^) 

with the projectors P^^'^ defined in Eq. (plj) . 

Taking the traces in Eq. ()96|1 for B Xuii^, we find that 

K = ][C[^^]\ (97) 



4 

71 



(n-q—n-q) 4 {n-q—n-qy 



^ 2{n-q — n-q) ' ^ {n-q — n-qY 



hi 



2{n-q — n-q) {n-q—n-q^ 
For S — ^ Xs7, where n-q = 0, the traces give 

ft; = ^(c!"-ic?'-cf)^ ft-, = o. ;4 = M(ci"-lc<'>-cfy. 

h; = -i (3C|" - 2C« - 2Cf - Cf ') (C|" - 2C« + cf ) , 

K = ^(cP^^?-d^r. (98) 

Here the /i°"'s depend on ln(/i/mf,) and in addition y = n-p/rrib through the 0{as) corrections 
to the Cj's (see Eqs. (ff2|) for results at /i = nib). 

The factors of ra-g = msx^ (for — > Xs7) and n-g = mB(l — Ijh)^ n-q = msil — Uh) 
(for B —>■ Xuiu) are purely kinematic and so can be directly replaced by these hadronic 
variables. Expanding to leading order about = 1 and uh = gives the results that 
should be used at LO in the SCET expansion. For notational convenience, we write 

= h^^ + /if ^ + . . . , hl = hT + /if ^ + . . . . (99) 

For B X^iu we have 

hr = \[c[^^]\ (100) 

zniByH f^EVH 
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while for B — > Xg'y we have 

hl^ = _l(3C'W-2C'f-2CP-C'f)(C'f)-2C'f +C'f), 

= ^{C?-lc?-Ci'r. (101) 



The resuhs in Eqs. pOO|) and ()101|) agree with Refs. |45|,|4a|. At NLO in the power expansion, 
it will be necessary to keep the next term from the expansion of the kinematic prefactors 
in Eqs. 1^ and For B XJu we find 

/i?'" = 0, (102) 



while for 5 we have 



<• - '"^'f '' (cM-lcW-C''y. A»'' = 0, (103) 





2 


-1) 




2 


-1) 


2 



5 2 2*^ 3/' 

In Ref. [sG], the LO triply differential rate for B X^fz/ was found to satisfy 
d?T^ / dxdydyo oc (x — ?/)(?/o — a^)- In terms of the this corresponds to the relations 
= Ah^^/yn and niBh^ = 2h^ /yu (where at LO we can set y = [l — uh){'^ — yu), 
yo = 2 — uh — yn)- Eqs. (|l(J(jp agree with this result at tree level in the hard functions, 
but give non-zero corrections to the first (/i^^) relation of order as{mb) from the Wilson 
coefficients Cg^"* and C^^^ . 

V. VANISHING TIME- ORDERED PRODUCTS AT 0{X) 

To work out the factorization beyond LO is now simply a matter of going to higher order 
in A in SCET. At order A the time-ordered products are 



f = -I jd^x e""-" T[j(°)t(a;) j{ia)(o) + j(^-)^(x) j(°)(0)] , (104) 
f.'^-./.*.e--^T[/«..W/-,0)W'-W.<«.(0)]. 
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We shall see that these time-ordered products give vanishing matrix elements for B Xs'y 
and B —>■ Xu^v. The result follows almost directly from chirality and the fact that the 
currents in Eqs. (|64p . (j65p and (j67|) all involve left-handed collinear quark fields. However, we 
must be careful to check that terms proportional to the chiral condensate do not contribute. 
We follow an approach similar to the previous section: by using the Fierz transformation 
in Eq. (jHT|) . we can collect the usoft and collinear fields into separate parts J ® and the 
matrix element then factorizes to give (OliyiO) ® {B^\S\B^) . These parts are still connected 
by indices and spacetime integrals, which are represented here by the ®. 

For T^" and T^^, we take the expressions for the currents from Eqs. ()64|) and ()65|) . The 
same usoft fields appear as in the leading-order T'^^\ and so the T-products involve soft 
matrix elements that are similar to the one in Eq. ()82|1 . namely 



{B,\T[n.{xrf^rXF^vf}nM] \B.) , 



(105) 



with similar matrix elements for T^'' but with T^!"'' 

J H' 

multiplied by the collinear matrix element 



T^l'^''. The first term in Eq. (jl05|) is 



J<^''^\0,x) = - (or {x^^V:^h{0)F^xiM] 

UJ \ L 



while the second term is multiplied by [j7'*-^"''(x, 0)]^. For T^* the analogous result is 

-1 



[X 



(106) 



(107) 



In both J'^^"'^ and the fact that we need an overall color singlet eliminates the pos- 

sibilities = 4, 5, 6. Owing to the presence of the a _L-index, rotational invariance also 
eliminates k = 1,2, leaving only = —'i^-y^/{2Nc). However, this term has the wrong 
chiral structure and also vanishes, i.e. 



xXXn ^ XuPrHIPlXu = , 



(108) 



so both T^*^ and T^^ give vanishing corrections. 

For T^^ we need to use the Fierz identity in Eq. (jHT|) twice to group together all usoft 
and collinear factors. This gives the usoft matrix element 



<S.|H.(a;)r5.?V,^ (z|>[,7?5) (l/)i^."rf 7^.(0) |i?.> 



(109) 



multiplied by the collinear matrix element of a four-quark operator, which we claim satisfies 



x';.{y)F^'xi 



[X] 



0) = 0, 



(110) 



To prove Eq. (jllOj) we begin by noting that color now allows T° terms as long as they occur 
in both F's, so either k, k' G {1, 2, 3} or k, k' G {4, 5, 6}. For either possibility the argument 
is the same, so for convenience we take the former case. Rotational invariance now demands 
that one of /c, /c' is equal to 3 and the other is 1 or 2. If /c' = 3, then the first pair of collinear 
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quark fields vanishes, just as in Eq. ()1U8|) . while if /c = 3, then the second pair of collinear 
quark fields vanishes. 

The results in Eqs. (|1U8|) and (jllUp rely on the underlying assumption that there is no 
structure in the vacuum that can flip the chirality. In QCD we know that this is not the 
case, since the chiral condensate (and instantons) take L ^ R. The above argument is 
valid because we have used chirality only at the //^ ~ QK scale where we are matching 
perturbatively, and not for /i^ ~ <^ QA, where it is badly broken. The same argument 
applies when perturbatively matching on to the weak Hamiltonian at /i ~ mw ^ A. 

VI. FACTORIZATION AT NEXT-TO-LEADING ORDER 

Since the order A contributions vanish, the first power corrections occur at order \^ = K/Q 
and will be referred to as NLO corrections. The NLO contributions to the decay rates have 
several sources, including the following: 

i) expansion of kinematic factors occurring in front of the Wi in the decay rates given by 
Eqs. and dng) , 

ii) expansion of the kinematic factors appearing in the hi at LO, i.e. in Eqs. (jU7j) and 

iii) expansion associated with the conversion from partonic to hadronic variables, given 
in Eq. , 

iv) expansion from higher-order operators contributing to the time-ordered products, 
given in Eq. pi8|) . 

Once the region of phase space where the SCET expansion is valid is properly defined, as in 
Sec. Ill C\ the corrections in i) and ii) are straightforward to compute. The corrections from 
i) do depend on the choice of how the xh charged-lepton variable is treated, for example 
whether we integrate over all of xh or instead look at an xh spectrum with a cut. Results 
for i) are given in Eqs. (j43|l . and (j^2|) . For ii) the required terms are derived by keeping 
one more term in the Taylor series in passing from Eqs. ()97|1 and fl98j) to Eqs. (jlOOj) and 
f)101|) . We give the results for i) and ii) in Sec. IVII Bl Note that these terms are already in 
terms of hadronic variables and are therefore unaffected by the conversion in iii) . 

The NLO terms from category iii) are discussed in Sec. IVI Al below. The contributions 
from category iv) require several sections. In Sec. IVI Bl we give a complete list of the time- 
ordered products arising at second order in SCET (category iv)), along with a summary of 
the jet and shape functions they generate. In Sec. IVI CI we carry out the tree- level matching 
calculations for these time-ordered products and define the jet functions at leading order in 
Us- In Sec. IVI Dl we give operator definitions for the shape functions. Detailed derivations of 
the NLO factorization theorems for these time-ordered products are presented in Sec. IVI El 
Finally, results for the computation of the traces that give the NLO hard coefficients are 
given in Sec. IVIFI 

Note that the final results for contributions from i), ii), iii) and iv) are summarized in 
Sec.lVlTl 
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A. Switching to hadronic variables at order 

At LO the factorization theorem in terms of partonic variables was given in Eq. ()94|) : 

j'~K—r+ 

W^^^ = h'^{n-p,mb,fi) / dk+ J^^\n-pk+,n) f^^\k++r+,i^). (Ill) 







In this expression, the variables n ■ p and n ■ r have a series expansion once we switch to 
hadronic variables. For the accuracy needed at NLO we have 

n ■ p = n ■ px — A , n ■ r = {A — n-px) — ^^^2 "^^"^^ ■ (112) 

Expanding Eq. (jlllj) gives the LO term in Eq. (j95j) plus the NLO terms 

j^,„,jNLO ^ _(^l_t^ fcO,p-_,„^, j^^. j,„,(p- , + ) (113) 

-A r''dk+ -^\h'^{n-p,mt)J^^\n-pk+)} _ /(°)(A;++A-p+) 







dn-p 



n-p=p^ 



(A1 + 3A2) 



2m 



B 



WXPx, m,) [J^'^ ® /(o)] {p\,Px) - A [{h'/j^'^)' ® /(°)] (p+ , p-^) 



Note that in taking the derivative with respect to n-p we differentiate only the Cj{n-p)^s in 
and not the prefactors depending on n-q. For later convenience, we switched to hadronic 
variables in the final line of Eq. ()113|) and defined 

[J^'^ ® m iPx.Px) = jfdk^ J''\Px k^) \^^l-litl±tzl^ (114) 



-5(A;+-p+)/W(fc++A-j9+)|, 



[(/i°^j(°))'®/(°)](p+,p^) = / dk+ ^\h'^f{n-p,m,) J^'^\n-pk-^) 







■n.-p=Px 



x/(o)(fc++A-p+). 



B. Time-Ordered Products at order 

To enumerate all the possible time-ordered product contributions at this order we consider 
all possible combinations of SCET currents and Lagrangians from Sec. IIIH 



where rti + . . . + rtj = 2 for NLO, i.e. (9(A^). It is useful to divide these time-ordered 
products into two categories, those that have a jet function that starts at tree level and 
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T-product 



Example Diagram 



Hard, Jet, and 
Shape Functions 



Usoft Operator 



f{2H) 




^ 




f{2a) 



f(2L) 




> 




x>--> V > — -6< 





;,3,4 o^(-2) .(4) 
''' •-'1,2 ^3,4 

;^3,4 (4) 



/i^,(x)(DT,±/i^)(0) 
(/i„Z?T,±)(x)/i„(0) 



/i,(x)(Z)x^±)(y)/i.(o) 



f (25) 




f(6) 
r(-4) „(6) 



•-'2-4 i/S-lO 



Kix)q{y)q{z)hy{0) 



TABLE IV: Time-ordered products that are of order = A/mh overaU, and that are non-zero 
at tree level. The power of is obtained by multiplying the powers from the jet functions 
by those from the shape functions / or g. We suppress color and Dirac structure in the usoft 
operators listed, which can be found in the text. The time-ordered product in the last row has 
not been considered in the literature and is enhanced relative to the other entries by a prefactor 
of AirasiExA.) ~ 5. 



those whose jet function starts at one-loop order only. To determine into which category a 
time-ordered product falls, we first note that the jet functions are vacuum-to- vacuum matrix 
elements, so all coUinear fields are contracted. Since there are no external _L-momenta, at 
tree level all coUinear lines have no ^-momentum and factors of V± all vanish. (This is also 
true beyond tree level for factors of V± that act on all coUinear fields in a J or J"!".) Thus, 
for example, the product j(i'')t j(i^) has a jet function that starts at one- loop order since 
we must contract both the coUinear quark and gluon lines. A second example consists of 
time-ordered products that involve a C^^^ insertion, in which neither the V± nor the A;^ in 
the can contribute at tree level. 

The time-ordered products that appear already at tree level will be most important 
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T-product Example Diagram ' ' Usoft Operator 

Shape Functions 




TABLE V: Time-ordered products that are of order A/mj,, but have jet functions that start at 
one- loop order. The last three rows introduce new shape functions that were not present at tree 
level. Vertices that are not labeled are from df). 
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T-product 



Example Diagram 



Hard, Jet, and 
Shape Functions 



Usoft Operator 



f{2Ga) 



f{2Lb) 




y^[2Ga] ^(4) -h,{x){D ^){y)K{Q) 



K{x)n-D{y)h^{0) 



j<{2La) 




■-(X) j(0)t 




K{x)D^{y)KiO) 



TABLE VI: Examples of non-abelian terms in time-ordered products that are of order A/m^, and 
have jet functions that start at one-loop order. 



phenomenologically, and include 



f(2^^) = -^ /rf4xAe-^^-^'T[j(°)t(a;)^£f (y) J(0)(0)] 

ji(2a) 

f{2L) 
f{2q) 



(116) 



K=a,b 



jd'xd'yd'ze-- T[j(°)t(^),^a)(^),^(i)(^)^(o)(o)] _ 



These time-ordered products give jet functions that are either the same as at lowest order, 
namely J^^\ or enhanced by two or four powers of A, namely J^^'^~^\ The corresponding 
shape functions are power suppressed relative to /*^°^ by two, four and six powers respectively. 

The most pertinent information is summarized in Table IIVI along with examples of 
Feynman diagrams and the forms of the non-perturbative usoft operators. The overall power 
of is obtained by simply adding the superscripts in the operator insertions. Dividing this 
overall power into that of the jet and usoft terms is slightly more involved since it depends on 
the individual operators. At tree level it is simply that every additional coUinear propagator 
enhances the corresponding jet function by A~^. Generating additional propagators requires 
inserting subleading operators with additional usoft fields that produce power-suppressed 
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shape functions. The power counting in SCETj restricts us to consider at most two 0{X) 
operator insertions, so this guarantees that the greatest enhancement from the jet function 
is 0(A~'^) and the set of possible terms is finite. 

From Table IIVI we see that for T^'^"^ and T^^"^ the coUinear fields give a jet function 
that is identical to J'^^^ defined in Eq. (jHH), but have A/Q-suppressed shape functions fl^2 
and /q^''. For T*^^'^-' we find an enhancement of Q/A in the jet function J'^~'^\ but shape 
functions that are further suppressed, namely f^^l^ which are down by A^/Q^.^^ Finally for 
2"(2i3) have the biggest enhancement, Q^/A^, through the jet function J'^~^\ and a shape 
function f^l involving a four-quark operator, which is down by h? /Q^. 

In Table nVl we also notice that it's only the J*^°^ currents that appear, plus the subleading 
currents J^^*^) and J^^^-* whose Wilson coefficients are related to the coefficients in J'^'^^ by 
reparameterization invariance (see Appendix El- This implies that the logarithms encoded 
in the running of the subleading hard functions h^~^ are the same as the logarithms that 
can be resummed in . (Here the term "logarithms" includes Sudakov double logarithms 
as well as the usual single logarithms.) This result depends on the fact that the Lagrangians 
that are inserted do not run. This implies that the logarithms resummed between the scales 
rnl and m^A are universal for these terms (besides the simple HQET running from cp{ii) in 
Cj^ , which is easily taken into account). Additional logarithms occur below the m^A scale 
and there is no reason to expect that they are universal. A computation of the anomalous 
dimensions of the corresponding soft operators would allow these additional logarithms to 
be resummed. Finally, for the time-ordered products that appear at one-loop order in the 
jet functions, even the logarithms between ml and m^A are not the same as at LO. They 
are, however, entirely determined by the running of the J^^-* currents that were calculated 
in Ref. 0. 

At one-loop order in the jet function, the remaining time-ordered products start con- 
tributing. They are 

f^^''^ = -I jd^xe-"-'' J2 T[j(i'')t(3,) j(i^)(o)], (117) 

K,=a,b, i=a,b 

f^'^^ = -il'd'xe-'^-^ J2 T[j(2^)t(^) j(o)(o) + j(o)t(^) j(2^)(0)], 
f^^^'^) = -z /"d^xAe-*""5^T[j(i'^)t(3.),£a)(^)j(o)(Q)^ j(o)t(^)^ 

K=a,b 

-i J d^x dS e-"-" T [ J(i'^)t (a;) (^) j(o) (g) + jWt (a;) (^) j(i-) (q)] , 

K=a,b 

f^^^") = -I jd^xd'ye-''-^T[J^'^^^ix)tCf^\y)J^'\0)] 
-% fd^xd^ye-"-^ T[j(o)t(a;),/:(2'')(^)j(o)(0)] , 



(2) 

Note that in the subleading Lagrangian C^^ the two _L derivatives appear at the same spacetime point, 
providing a simple explanation for why this occurred to all orders in the twist expansion in Ref. 5l|. 
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-t J d^x dS d^z e-^'^-" T [ J(°)t (3.) ,£(1) j(o) ^0)] 

jd'xd'yd'ze-- T[jWt(^) j(o)(o)] , 

f^^^'^) = y"d^a;Ae~*""T[jWt(a;),42«)(^)j{0)(0)]. 

Information about how these time-ordered products contribute to the factorization theorems 
at 0{\^) is summarized in Tables |3 and IVIl Note that T'^^''^ and T*^^'^-' match on to sub- 
leading jet functions and the leading-order shape function, and so these 0{as) corrections 
can be calculated without introducing new non-perturbative information. The remaining T- 
products give new subleading shape functions that were not present at tree level. Together 
the terms in Eqs. ()116|) and pi7p provide the complete set of time-ordered products that 
contribute at this order in A and at any order in a^. 

Adding up the contributions from the various subleading time-ordered products, we find 
that 



= ^^^^ rdk^ J^'\n.pk\^,)f^\k^ + r\^,) (11 
^rrib Jo 

+ V^^^ rdk+ J^'\n.pk+,^) jf\k+ + r^,f,) 

r=l Jo 

+ ^^i^ [dktdk^ j}^?{n.pkl,f,)n'\kl + r\f,) 

r=3 

+ ^^SM [dkfdktdkt Ji-'\n-pkl„f^) e\kl, +r\^^) 
+ ^^ISl^ rdk+ J^''\n-pk+,fi)g^^\k+ + r+,fi) 

ITT'b Jo 

+ ^^i^ [dkfdkt Jt!\n.pkl,^,)gi'\kl + r\^,) 

+ y"^£iM jdktdktdkt Jt'\n-pk^,,^i)gf\k^,+r+,^i) 

+ ^^i^ fdkfdktdkt [Jt"\n■pk|.^^)gi'\kf + r\^^) 

+Jt'\n-P k}, ^^) gl%{k} + r+, /i)] 
+ J2 dz.dzA ^^ll^^l^ / dk^ J^\z,,z,,p-~^k^) f^'\k^ + K-p\) 
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^t^4 Jo 

10 „ ,[2c]m+8/ ^ p+ 



m=5 

[2LL]f^U) 



ISA J ! 



^2) (—2) (—2) 

where j = 1,2 and j' = 1,2,3 and J'i±2 = i <^2 for r = 3,4 respectively. Recall 

that f = s for B — > and f = u for 5 — X^^z/. If we work at tree level in the jet 
functions, then only the first four lines of Eq. (|118p are required (with shape functions /o-e)- 
In our notation, the fi denote shape functions that appear from tree-level matching, while 
the Qi denote shape functions that show up only when one-loop corrections are considered 
for the jet functions. 



C. Tree-Level Matching Calculations 

In this section, we show how computing the four tree-level diagrams in Table IIVI im- 
mediately allows us to obtain the tree-level jet functions, properly convoluted with the 
non-perturbative shape functions that show up at NLO order. Note that no summation of 
operators is necessary. An operator-based derivation of the factorization theorems at all 
orders in the jet functions is given later, in Sec. IVI El 

First consider the computation of Fig. Elusing the LO jet function from the time-ordered 
product in Eq. (f75jl . with currents as in Eq. (jHUj) . Working entirely in momentum space and 
contracting the coUinear quark fields in the diagram gives 

{-^mnT^'\^ ^_ . , (119) 

^ I [n-pK+ + It) 

where n-p is the large momentum flowing through the coUinear quark propagator. Taking 
(— I/tt) times the imaginary part gives 

jr(0)(A;+) = 5(A;+) + 0(a,) , (120) 

for the tree-level jet function, which is equal to 5(£^ — r+) by momentum conservation since 

£+ = A;+ + r+. (121) 

Thus, from Eq. ()119p we can directly write the convolution of jet and soft factors as 

Tr[p4°^|r(°)] j dt{j^\t-r^)] {KY5{t-in-d)Y%,], (122) 
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FIG. 4: Momentum routing for a) T^^^^ and b) T^'^'^\ In a), the operator with one gluon comes 
from [i-^Ms]- 



where we have used Ti^ = Y'^hy. Performing the integral at tree level gives the well-known 
fact that the decay spectrum is determined by the shape function evaluated at r+. Note 
that in SCET we needed to compute only one diagram rather than resumming an infinite 
sum of operators. 

At NLO the tree-level diagrams are shown in Table IIVI The graphs for T*-^^-* and T*^^"^ 
involve the same calculation as in Eq. ()119|) but just leave different soft operators. 

For T^"^^^ the momentum routing is shown in Fig. Contracting the collinear quark 
propagators in the graph gives 



n-p 



-Tr 



2 [n-pk2 +ie 



J {k++ie){k^+ie) 



n-p 2 2 {n-pkf +ie) 



(123) 



where momentum conservation sets kf = if — r"*" and fc^ = — r"*". Taking (— I/tt) times 
the imaginary part gives the tree-level jet function 



— 



(124) 



and the convolutions 



n-p 



-Tr 



x\KY6{£t - tn-d)YUD^^''tD^f^Y6{it - in-d)Y^K 



(125) 



Finally, we consider the NLO time-ordered product T^^'^) with the momentum routing in 
Fig. HId. This graph gives 



2 {n-p fcg +ie 



7lT^^„.(£^-4) 



n-pk2+ie 



(126) 



9 



n-p {ki +ie){k2 +ie){kf +ie) . 



2 (n-p k{+ie 

M 2 
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where kf = ii— r+, fc^ = £^ — r"*", and = — r"*". Taking (— l/vr) times the imaginary 
part gives the tree-level jet function 
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6{kt 



+ 



+ 



7r'6{kt)6{kt)6{kt 



[{ktm) (ktm) (ktm) 

Writing out the convolutions for these four-quark operators, we find 



n-p 



(127) 



:i28) 



D. NLO Soft Operators and Shape Functions 

In this section, we define the shape functions that can appear in the NLO factorization 
theorem. At tree level, many of these functions have already been defined in Ref. (we 
use a sequential enumeration of these functions, and translate to the notation in Ref. (51| 
in Table IVTT]) . We also define the additional shape functions with four-quark operators that 
are not given in Ref. [H^ , and later in section IVI El the shape functions that appear only 
beyond tree level in the subleading jet functions. 

We begin by enumerating the usoft operators. For the time-ordered product of the LO 
shape-function operator and the subleading HQET Lagrangian we define 



Of\t) 



'dx 



e 2" 



-1+ 



in 



d'yT[K{x)Y{x,0)K{0) tOniy)]. 



(129) 



With igG^g^ = [iD:^^,iD:^g], the remaining ultrasoft operators that play an important role 
include 



Of(£^ 



dx 



Stt 



e 2 



n,ix){Wi^M{0) + (h.c, x^O) 



1, 



2 

dx 



K{tDi,6{e+-in-D^,)}K 



(130) 



in 



e 2 



e+3 



n,{3:hii5{^in,m + (h.c, x^o) 



1 ^ 



K{tDi,6{£+-in-D^,)}^^-f,K 



dx 
87 



e 2 



e+3 



-iTH^{x){iDf^^n^){0) + (h.c, x^O) 



K[iD^^^,5{t-in-D^,)\K 



dx 
StT 



e 2 



t+x- 



-tTn,{3:hl^,{Win,m + (h.c, x^O) 
-m-L'us)l7T75^« , 
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= hJ{i^-in-D^s)9G"'^^6{if~in-D^s)lll5hv , 
and the four-quark operators 

OfiiUs) = f"\fj"~ e-'^^'t'^-^^^^ (131) 
= I {M(^3+-m-Z}„,)7^PiTV} Siit-in-d) {gV^L5(^^-m-D„,)T^/i4 , 

= ^ [/i.5(£+-2n-Z}„,)7^Pig"] [gV^L5(^|-«^-^«s)M • 

Here V'us = {'^'4)/^'^us and the flavor of the quarks in these operators is s for B Xs'j 
and u for P X^^v. Note that a minimal basis of Dirac structures for the bihnears is 
/i^{l, 7^75}/ii,. In the second hne for each operator we have used our freedom to integrate 
by parts since only the forward part of these usoft operators is required. The o per ators 
in Eq. ()130|) were obtained by matching with tree-level quark propagators in Ref. j5l| and 
we agree with these. The four-quark operator O"^ appears when gluon propagators are 
included at tree level, as we discuss further below. The operator J difdi^OQ^ occurs from 
the disconnected annihilation contribution, as shown in Ref. jHll. However, since this con- 
tribution is of order (A/m^)^ in the endpoint region, we do not include it in our analysis 
here. 

When 0{as) corrections are included in the sub leading factorization theorems, then in 
general we require additional usoft operators. For example, at 0{as) there are two additional 
soft operators generated by T^'^^\ and six additional operators from T^'^'^\ Definitions for 
these operators can be found in Sec. IVI El 

The subleading shape functions of the operators in Eq. ()130p are defined by the matrix 



It seems that Ref. |^ has a typographical error in the phase factor of O^'^ and and that it should 
read e*'^'^i~'^^^*ie~"^i*=' , in which case the variables in that paper are related to ours by uJi = —ii. 
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elements 



{B,\Oo{t)\R:) = fli'\t), (132) 
V n-v/ 

(5,|Pf = 6f /f (£+), 

{B^\Ofin\B^) = gfft\it,^t), 



'1 ' ^2 



where we take vj_ = and = e^^'^'^v„nT-/{n-v). For the four-quark operators we need 

n^n^{B,\Of{£l2,3)\B.) = fti^JtJt) , (133) 

i9a,-^eip){B.\0f{il,,s)\B.) = fti^Jt^K) ■ 

Recall that we use the notation /j to indicate shape functions that appear from tree-level 
matching and Qi to indicate shape functions that show up only when higher-order perturba- 
tive corrections for the jet functions are considered. Definitions for all the Qi functions are 
given with the derivation of the factorization theorems in Sec. IVl El For the /j, we give the 
translation to the notation in Ref. .SlJ in Table I VI II 

In Ref. jsj], it was shown that the equations of motion for hy imply that 

/f)(.;) = 2^/^(0.) , (134) 

which reduces the number of NLO unknowns, and that the B matrix element of the operator 
O2/3 vanishes, so we shall not need it for our results. The same is true of two additional 
operators, O^^^ and ^3"^, which are the analogs of P^'^ but with anticommutators. Ref. [HH 
also removes from the basis by noting that all the moments of the corresponding shape 
function vanish. However, we keep this operator since it is unclear whether the moments 
are sufficient to define the function completely beyond tree level in the jet function. (We do 
show, however, that this operator is not matched on to at tree level in the hard function.) 

E. Factorization Calculations at 0{X'^) 

We now derive the factorization theorems for the subleading (9(A^) terms. We show that 
it is exactly the soft functions defined in the previous section that appear in the decay rates, 
and give operator definitions for the jet functions that can be used beyond tree level. 

1. Calculation of f^"^^^ 

For T^'^^^ the factorization is identical to that at LO in Sec. llVt except that the final 
usoft matrix element is the time-ordered product 

j'^e-i'--j,'y{S,\T \k.mmM .te. (136) 
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Our Notation Notation in Ref. 



/ I*- ) 




JQ I*- ) 


f( f+) 






/f (£+) 










h2{-it,-4) 











TABLE VII: Relation of our /, functions to the notation in Ref. 



The factorization theorem for T^'^^^ is therefore 

lyf/^) _ h^^(n-p,mb,fi) r dt J^^\n-p{i++mb-n-q), fi) ^^^^^^ . (136) 

Jn-q-nib Znih 

As noted in Ref. the function /q^"* can simply be absorbed into the LO f^^\ since at 
this order they always appear together in the combination 

/^°nn + ^/o^'^(n- (137) 

2. Calculation of f 

For the remaining time-ordered products, the factorization is more complicated. How- 
ever, some aspects of the LO analysis in Sec. IIVI remain the same at subleading order. In 
particular, we still have the exp(— ir ■ x) phase factor in Eq. (f?^ . a 6{uj' — n-p) appears in 
the separation of hard Wilson coefficients, as it does in Eq. ()79j) . and the vacuum matrix 
element of coUinear fields still gives a 5{uj — uj'). For convenience, we integrate over uo and 
oj' and remove these delta functions right from the start when considering T*^^"\ 2"(2i)^ g^^^j 
T(2g)_ 

For f the analog of Eq. ^ is 

Tt""^ = -^E E CAn-p)Af''\n.p) jd^xe^- (138) 

3,3' K.=a,b 

XT [Jfj{n.p,x) jf:\n.p,Q) + jf;^\n.p,x) 4°i(n-p,0)] , 
where aJ''^\uj,ij) are the Wilson coefficients for the currents, as in Eq. (jUHj) . The 
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products of currents in square brackets are 



.(0)t j{2a) j{2a)] j{0) 



1 



2mh 



(2a) ■'^Ta ■ 



(0). 

r 



jmj{2b) _^ j(26)t j(0) 



(26)„ 



+ -i [(K,(-0^;^r)Tj.f (x) [wr;?^^.] (0) . (139) 



(0). 



Next we insert identity matrices to Fierz transform using Eq. f)81|) and take the (5^ 
matrix element. These two terms then become 



(0)t j(2a)^ji2a)]ji0)\^ _ ("l) 



B^\n.{x)Tf Fj:Tf:\i'Dl':nm\B.) (140) 



-(2a) f-'^Ta 



[B,\{nM:)ix)r%^''F^rpnM\B.)]{o\xnAo)F^^^^^^ 



I jm j{2b) ^ j{2b)] 



1 



I {x)vfF^Tf^^ 0^^n.m I B,) 



An-p 

-(2^) Fnr(O) 



\B,)/2 



-{B,\{n.ivi:){x)i:)jF^TfnM |5,>|(o|xn,.(o)F,^xnx(^)|o> , 

where it is understood that we have u = oj' = n-p, and we have used g_L = in order to 
integrate the iD-^ derivatives by parts in the J*^^^^ terms. From Eq. ()140|) we see that the 
colhnear vacuum matrix elements are the same as in the LO case. Hence they are non-zero 
only for k = 1 and are given by the LO jet function J'!^\k^) from Eq. (j83|) . 

Using the trace formula in Eq. (jHUj) . we can simplify the spin structures in the usoft 
matrix elements in Eq. ()140p . where we can now set F" = ?i/2. Here both the and 
terms in the matrices will give non-zero contributions, so A;' = 1, 2 and we have 



9a/3_ry{2a) .{2a) da/S ry{2b) .{2b) 



[B.,\V'H,{x)Pi:,{iVt'Hm\B.) 



-Tr^P,^ 



9aP_-Y{2a) . (2a) _ 9a£_-y{2b) {2b) 



^ Tf]{B,\T {n^ivD {x)p^,nM\B.) 



= -\i:v{x^;=']m{o{ + tO',)m + h:T{x'^='] 

-^Tr{X,^;=2}(5.|(Pf^ + .Pf )|P.) + ^Tr{X,t'}(5.|(-Pf^ + .Pf )|P„) 



(141) 



where 



X' 



X 



I T ± 
H^{^) t \ 9aP_y{2a) .{2a) _ 9al3_,^{2b) .{2b) 



p^F' — 

^k' J- 2 



pH 

^k' 



9al3 Y(2a) (2a) _ 9ap ^(26) {2b) 



2 



(142) 



Note that the second trace is the complex conjugate of the first with (/x ^ z/), since 
loiF^Ylo = and Tr[X] = Tr[Xt] = Tr[X*] = (TrX)*. These usoft matrix elements give 
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the shape functions /i, go, and /2 in Eq. ()132p at any matching order in as- The index P is 
transverse to v, as follows from Eq. ()14H1 with v± = 0. Combining Eq. 1)14111 with the delta 
functions and prefactors from the jet function and the (— i) from T^"^"-^ gives 



(143) 



x{2)/,+ I ^+ 



r 



We include these traces in the definition of the hard coefficients at subleading order, as was 
done in Eq. ^^i . defining 

^(0) 1^ / rf 



h 



[2a] 1 



1 ^c, 4" T.{pj;»i|( -^)t<-' + (h.c. ,».)]p, 



i ' 



(144) 



J,3 



-(0) i li€j^r^{2a).{2a) _ (2b) .i2b) 

. 4 ^ 2n-p ^""^ ^ . 



+ (h.c. ^ z/)| 



,[2a]00 ^V^r^ / D TDr(°)^ 



e^^(2a) .(2a) _ £l_^^(26) . (26) 

4 ^' 2n-p ^'^^ ^ 



+(h.c. /i^z/)}i^^r 



For the Dirac structures T*^^"^^ and T*^^^^ that are present at tree level is zero, so g'^^ 

comes in with as{mb) suppression. We were unable to prove our suspicion that g^ would 
not appear beyond this level. 

Using Eq. ()144j) with Eq. ()143)1 and taking (— l/vr) times the imaginary part gives the 
NLO factorization theorem for Wp""^: 



W} 



rUb 



(145) 



+ 



+ 



rrib 



'dk+ J^°\n-pk+,fi) fP{k++r+,fi) 



hf''^°°{n-p,mb,fi) 



nib 



Explicit results for the coefficients hf""^ are summarized in Sec. IVl Fl 



3. Calculation of r(2i) 

For T^"^^^ the heavy-to-light currents are identical to those in T^^\ and so the same Wilson 
coefficients appear as in Eq. ^7^ . i.e. 

= -^5^Q'(^-P)Q(^-P) Jd'x jd'ye--T[jP;i{n.p,x)Cf^\y)jf}{n-p,0)\ . 
hi' 

(146) 
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The difference is the extra Lagrangian insertion with mixed colhnear and usoft fields, namely 
jf\x)Cf;\y)jf{0) (147) 



k,k'=l 



X [xn{y)F^'Vn,Ax)] [xLio)Fi^''pXniy)] , 

where in the second line we have used analogs of the Fierz formula in Eq. ()81|) that are 
appropriate with left- and right-handed projectors, 



Pr. 



2Nr 



k=l 



R 



1 = o E^^^'^^fc""" 



k'=l 



4 



ill 



(148) 



The vacuum matrix element of the colhnear four-quark operator in Eq. p47|l must be 
a color singlet, implying that either k,k' G {1,2}, so that both P^'s have no T"'s, or 
k, k' G {3,4} and the color structure gets reduced when we take the matrix element, using 
T°- ^ = ® T"^ 6"'^/{N^ — l). Rotational invariance then implies that the matrix 
element gives 6k,k' so we do the sum over k' and set k = k'. This leaves four terms in 

m « m - 4, « [^pj + <^i±^ [Mpj « [Mp,] ,149) 



-2Nr 



-2Nr 



^[^PlT^]^[^PrT^] 



When we use Eq. (j86|) to reduce the Dirac structure between the heavy quark fields in 
Eq. ()147|) . and take into account the (5"^'s from Eq. ()149|) . the usoft operators become 



J2 Tr P.^r5?V,^"7l7lfi^rM°^ n^{x)F,^{^Vi^^V^:){y)nM 



(150) 



k"=i 

2 



I E Trter5.?V,^-7l7l^i^rrf 1 [^.(x)P,^T^(^2^iriPi-;)(y)T^H.(0)' 



k"=i 



for G {1,2} and k G {3,4}, respectively. The spin-trace that occurs in both terms can be 
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simplified with 7^7^^ = {gY ^^^Y'l^)^ ^llll = 4{9±^ ~^^'i^'l5) 



Tr 



Tr 



k"^j' 2^ j 



(151) 



The — ie" ) terms contract with the ida/s+'^^a/s)/'^'^ Eq. (|149|) to give a 2. 
From Eq. p49p there are four jet functions: 







(2vr) 



2 _iri,+ 
e 2 



(152) 



2(0 



2NcV 



Xn{y) 



uj J [2tiY ■> 



(iV2-l) 



XnmPLT^nA^)] [xt.mPLT^^Xn{y) 



5\x^)5{x+)5\y^)5{y' 



rdk+dk+ ^ll,+x-+(k+-k+)y-] M~2)^,+. 



(iV2-l) 



XniynilPLT''Xn,.'i^)\ [Xn,.m^tPRT''pXniy) 

dk-^ dk<2 rz-."i" 



^^^"^ '"'U\x^)6{x^)6\y^)6{y^ 



UJ 



{2n) 



e 2 



[k+x-+{kt-kt)y-] jl^^\k+) 



(—2) 

At tree level, only J'^^ is non-zero. It has the value 



jL^\n-ki,n-k2) 



[n-ki + ie] [n-k2 + ie] 



(153) 



At one-loop order, J/^'"^^ will become non-zero, while and J'^~'^^ are non-zero only 

at two- loop order. The expansion here is in as at the jet scale, /x^ ~ rux- We denote the 
imaginary parts by 



-1 



TT 



r(-2) 
'juj 



(154) 



Eqs. p49|) and also show that if we work at any order in ag there are four shape 
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functions generated by T^'^^\ i.e. 

\ n-v/ 



n-v 



1 ? ^2 ; 5 



{9l + ^el) |^^e-t[-^^+^-^^] {B^mx)Fj:,,T^{^Vi:^V^:){y)T^^^^^ 

= 4"i gMAt) - Sk"2 {vv-^)9^i^tJt) ■ (155) 

Here the index t] is from F^ = Pjj'j'^'j^Py. At tree level in the jet function, only the first of 
these four combinations occurs. The functions fs^^ were defined in Eq. fll32j) . The definitions 
of 5f3,4 are 

{B,\P^^{t)\B^) = -ef (£+,4), (156) 

where 

Of {it, it) = j%^e-¥^--e-'^^'t-^t^y-^.{x)T^^^^^ , 
P,f{it,it) = ^ f-^^e-¥^^-e-m-4)y-ln.^^)T^[^^^^^^ 

(157) 

The hard coefficients can now be defined to be 



h 



1-13 _ -^j:c,C,Tr{^ptp^r}pr--^f^, (158) 



n-p — ' I z z ■'J n-p 



n-p ^-^ n-v I / ■'2 ■' J 

3,3 
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Combining equations, we get the final result for the T^^^^ contributions: 



^[2L]3/- X 



,[2L]4/_ N p 

+ ' ^2^' ldktdkt{[jt'\n-pkl)-jt'\n-pkl)]Uikl+r^) 

+ [jt^\n-pk+) - Jt^\n-pk+)] (74(A;++r+)} . (159) 

Explicit results for the hard coefficients /if ^' are summarized in Sec. IVI Fl 



4. Calculation of T'^'^'i^ 

For T^^'^) the heavy-to-light currents are also identical to those in Eq. (f7S|). but now we 
have two insertions of the Lagrangian £^^^ , i.e. 



f{2g) 



\Y.C,,{n.p)C,{n.p) jd'x JdSe--T[jfJ{n.p,x)Cf^{y)Cf^{z) jj?(n-p,0)] . 

(160) 



Owing to the hermitian-conjugate terms in each jC^^g , there are two ways to generate the 
four-quark operator shown in Table IIVI where the positions in the figure are labeled either 
— z — y — xotO — y — z — X. This average over z ^ y cancels the 1/2 in the definition 
of T^'^'^\ Fierz transforming the product of operators with Eq. ()148|1 in order to group the 
collinear and usoft fields, we have 



jf\x)qj{y)C\\>{z)J^'>iO) (161) 
= I E f(xn^^?^l)(l/)^r'4.'(a:)lx^,jO)F,^^(-^^rxXn)(^) 



'(1)/ 



'(1)/ 



k,k'=l 



X 



The vacuum matrix element of the collinear operator must be a color singlet, so either 
k,k' e {1,2} or k,k' G {3,4}, and then rotational invariance implies that k = k'. This 
leaves four non-trivial collinear operators, which are given by 



(162) 



with the decomposition of F]^^ ® F^"^ given by exactly Eq. ()149j) . For the usoft operators 



we always have 



so this projects on to two components of the light quark fields 
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given by ipus with 



l.n 

us 



0. The Dirac structure can therefore be simphfied by noting that 

^PtPL,Rrn, = Tr[^P,.,^rP,]Cs^L,/?|H. + Tr[7^^rP„]C.i'L,ij7l^., 

H,rP^,iC, = TT[TPR,L^P,]n,^PR,L^t + Tr[r^7.^n]^.7li^fl,LCs , (163) 

and using rotational invariance of the usoft matrix element to determine that 7^ and 7^ 
terms are restricted to appearing together. Thus [Hv'^j' Fk""^us] [i^usF^^T^J'^H^'j together 
with the (5'Q/3+^e;^^)/2 factor from Eq. p49|) gives 



Jkl 



Tr 



[7Y.^p«C] [As'^PRn.] 



+Tr 



r;?^^.^7in]Tr[7lfp.rfn 



-s, 



fc2" 



r — 2 



[n^irPni^l] [^li^PRn.] 



Tr 



TPfr.PnP.]Tr[^llPnrrP. 



[^v^Pl^U [i^l^PL-H.] 



+Tr 



rPPL^7lllP.\Tr[^llllPnTfP.\ [^.T.^^lC] [C.T.^^l^.] 

} +44{-- 



Jk3 



(164) 



where the structures for the 6k3 (5fc4) term are identical to those of 5^1 {Sk2) except with 
extra T^CST"^ factors in the operators. Simplifying the Dirac structures and using rotational 
invariance of the (Pt,| ■ ■ • \Bjj) matrix element gives 



Ski < Tr 



r(o)#p p 



Tr 



^PnTfP. 



[n.^PRrus] [rus^PRn.] 



+Tr 



rf^7lP.P.]Tr[7lfp.rfp. 



-S,2i Tr 



' 2 

J- j' -ttIlPrPv 



Tr 



'-^ I p p{o) p ■ 



+Tr 



rJ.?Vp^p,lTr[^p^rf p. 



[n^^^PL^i^t] [i'tilPL'H.] 



+5, 



k3 



} + 5m{ 



(165) 



From Eq. p62j) the jet functions are 



4 



AiNs fdkfdk:^dk 



UJ 



2N, 



3 „-l[fc3+x-+(fc+-fc3+)y-+(fc+-fc+)^-] J^-'^\k+)^ 
illPL 



g 2 ^"-3 



(a;)x„,Jo)- 



4ziV5 fdkfdk:^dk^ 



UJ 



i2ny 



^l[k+x-+{kt-kt)y-+(k+-k+)z-] J^~^\k+) , 
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AiNs fdkfdk+dk^ 



(2vr) 



AiNs f dkf dk2 dk^ „iri.+o 



00 



(2vr) 



l[k+x-+{kt-kt)y-Hkt'k+)z-] J^-^\k+) 



1 > ' 



(166) 



where the common prefactor is A^^ = 5{u — u')5'^{xi_)5{x^)5'^{y_[_)5{y^)5'^{zi_)5{z^) and the 
1/V factors in the operators do not act outside the square brackets. At tree level, only 
Ji^'^\k^) is non-zero. It has the value 



jL'^\n-ki,n-k2,n-k3) 



[n-ki + ie][n-k2 + ie][n-k3 + ie] 
As before, we denote the imaginary parts by 



(167) 



TX 



r(-4) 



(168) 



From Eq. ()165|) there are eight shape functions with different Dirac and color structures, 
namely 



'dx dy dz 



.g-i{4.-+(4-4)j/-+(^^-4)--}(5,|?^,(x)^p«^^^(^)^^^(5)^p«7i,(0)|s,) 



1287r3 



{B,\n,{x) iPRT^^:M^Uz) iPRT^nM\B.) 



1 ? -^2 ? -^2 ; 5 
'dx~dy~dz~ _i[ ... j 
1287r3 ^ 



dx-dy-dz- ^.^ - l(^.l^.(5)7lPi?T^C(y)Cs(5)7^PijT^^.(0)|5.) 



1287r3 
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= /6(£+,4,4), (169) 

where ipu^ = (^^) / 4'?/'us and the elhpses denote exactly the same exponent as in the preceding 
expressions. Only the last two of these shape functions show up when we work at tree level 
in the jet functions. 

We define the hard coefficients to be 



^1* - iSEC^/CiT.{pj5?'fp,}T.{p„fp«rf}p, 



(0) \ p^,u 

i ' 



,,,8 =^-^J2 Cf C, ^Pjf \^-P,] T.{p.ll\PnVf)pr , (170) 
and the final result from the T^'^'^^ contributions is then 

J2g]7 



h r ( ^ 

Vrf'^) = ^ j dktdktdkt[jt^\n-pkj)g^{kj+r-^) + jt^\n-pk-^)gj^^^ 

dktdk+dk+{^jt^\n-pkpgwikf+r+) + jt^\n-pk+) gs{k+ +r+)^ 
+ J dktdk+dk+\^jt^\n-p k+)g5{k++r+) + ji~^\n-p k+) h{k++r+)j 

dktdk+dk^^^jt^\n-pk^)geikl+r+) + jt^\n-pk^) h{k+ +r+)^ . 



J2g]5 



+ 

nib 



(171) 

Explicit results for the hard coefficients /^p''^^'^ are summarized in Sec. IVI Fl 



5. Calculation of f (^^^ and f^^^^) 



For later convenience we set 



j duj^dcu2 Bf\uoi,uo2) Jf\u:i,u:2) , (172) 



where 



Jf\^i.^2) = ^Xn,.A^9B:^U , (173) 
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and Bj^^-^ {uJi,uJ2) = —Bj''"''^{uJi+uj2) mb/{uJi+uj2) for j = 1-3 when f = u, and for j = 1-4 



when / = s, and Bj {uji,uj2) = B^^"'\ul)i,ul)2) I {n-v) for the remaining j's. 

Both T^^^^ and T^^^^ have jet functions that vanish at tree level but are non-zero at 
one-loop order. The steps are the same as in previous sections. Here we have 

fi^b) ^ _^ ^ jdujiduo2du'j,duji i?]^^(u;i^2)-B]P(co'3,4)5(cj3+cj4-n-p) 
3,3' 



(174) 



To Fierz transform we can use 



4 



Pl®Pr= {^Pr) ® {^Pl) + (|PijT») ® i^P.T'^) 



(175) 



and we can drop the second term, since color implies that it does not contribute in the 
matrix elements. Thus, 

-1 



T{l)t t(1) 



2ml 



Tr 



Pv ^(1) i p ^(i)' 



;i76) 



X] 



where we have used the fact that the second term in the reduction in Eq. ()8(j|l gives a 
vanishing matrix element. Between B mesons, the soft operator in Eq. ()176|) gives the 
leading-order shape function f'^'^\ The matrix element of the coUinear operator gives the jet 
functions 



dk+ 
27 



e 2 



(0| [(w.^^?^i:J(0)^i^L(-^^7Si::Xn,.3)(^)] |0> (177) 
= {-2i){uji + UJ2) 5{ijJi+uj2-oJ^-uJi}5'^{x^)5{x'^) 

and their imaginary parts are denoted by J'm\ijjjk'^) = (— l/7r)Im jT^i^j (fc"^)- Noting the 
delta functions in Eqs. ()174|1 and ()177|) . we define the hard coefficients to be 



/.P(.„.„n.p) = J2B!^^\us,.)BfM9T'^Tr 

< ' ■> ■' m, 



3,3 



nib 



hf^'\z„Z2,n-p) = $:5«(..3,4)5f(..i,).6r'^Tr 



3,3 



TTlh 



2 j't^'fj- 2 joi'^ 



pr, (178) 



i 1 



where coi = Ziu-p, 002 = (1— Zi)n-p, 003 = Z2n-p, = {l—Z2)n-p. This gives the factorization 
theorem 

hr%uZ2,n.p) 



(2b) 



dzidz2 



rrib 



dk+ jr (^1, Z2, p^ fc+) / W (A;+ + A-P+ ) (179) 



hr^'\z,,Z2,n.p) 



+ dzidz2 



nib 



Px 



dk+ Ji'\zuz2,p],k+) f'^'\k+ + A-p+). 
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For T(2'^) we have 



(0), 



(180) 



where we have n = 1 for 



c and n = 1,2 for 



,j ^~.„^ C^^f^T^.^^^?-^^, so that C^^^^ contains the product of colhnear fields in each current. 
Fierz transforming T^^'^) with Eq. ()175|) (again we need the first term only) and then reducing 
the Dirac structure using Eq. (jHB|) gives 



d,e, f. For convenience we set 



jmjm j{2^)t j{o) 



(181) 



tr 



+ — tr 
4 



12^2 



^Tr^p«rl?) + h.c. 



Tr^Pi?r;.?)-h.c. 



n,{x)nM\ [xn,.(o)^cSl,(x) + (h.c, x^o) 
n^{x)nM] [xn,.(o)^Cl(^) - (h.c, x - o)' 



The usoft operator here gives the leading-order shape function f^^\ Next we define the jet 
functions 



-l){0\XnU^f^{[^t9n-B,]xn) ^{x) ± (h.c, X ^ 0)|0> 



-2i 



6{iJ—LJi)6'^{x_i)6{x' 



— (0\xnU^)^^[{'i'9n-Bc)-^2Xn,uji]{x) ± (h.c, x ^ 0)|0) 



-2i 



5{uj — UJi—UJ2)5'^{Xi_)5{x'^ 



-1 



(o|w(o) 



rUb 
iPL 1 



dk^ ±k+x- q-m 



2tx 



[z(?Sil^(zP,"^Xn)^J (a:) ± (h.c, x ^ 0)|0> 



—2i 

b{uj — ijJ\ — ijJ'2)b'^{xAb{x 

nib 



2tt 



7t<;i,2,8wi,2 



— (0|xn,.(0)^-[(z^7S,^"'^P,"J.2Xn,.J(x) ± (h.c, X - 0)|0> 
nib 



-2z 
nib 



6{uj — iUi — iU2)S {x±)5{x~^ 



27r 



9i^i o,10a;i 2 



(182) 
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and their imaginary parts Jm\oJjk'^) — (— l/7r)Im^t^^. (A;+), as well as the hard coefficients 

- ^^Y.^r{n-P)4'\n-P)^\^^^^^^ (183) 



2n-p ^ 



where uji — zin-p and 002 — {l — zi)n-p. 
Putting all the pieces together we have 



(2c) 



E 

m=3,4 
10 



,[2c]m+8/- X 



m=5 



hf^]^+^(^zi,n-p) 
nib 



/ dk+ji'\zuPxk+)f'^'\k++A-pt). 
Jo 



(184) 



^. Calculation of the remaining T^^^ 's 

The remaining T*^^^'s for which we have not yet defined hard, jet and soft functions 
include T^'^^°'\ T^'^^^\ 2^(2LL) rp{2Ga) ^ These time-ordered products all have jet functions 
that start at one- loop order and, unlike T^^'''^'^^ they also induce new shape functions, gj. 
Factorization formulae can be derived for these contributions by following similar steps to 
the previous cases. Rather than going through this exercise, we instead list some of the soft 
operators and shape functions that would be required. 

For r(2-^") we have 

{Pv\Ps,9{^l,2)\Pv) = ^± 5'll42(-^l') -^2") ) 

while for T^^^^) we find 



(185) 



(186) 
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For both T^'^^"-^ and T^'^^^\ the induced jet functions are of order A°. 
For T(2i^) the necessary soft operators are 

<A(^t2,3) = ^ f'^ulJ'' ^''-^ ■■■^[^.(5^)^^^i-r(y)^^^„^/(5)7r75H.(0)] , (187) 

together with operators Oii_i^{if2 3) -Pii-i5a(^^2 3) ^^at have different color contractions 
between the gauge-invariant objects Hv, iV^"', iV^'^, and Hv- The shape functions are 

{B.n'.i,xin,2,3)\B.) = -ef{v,-^)g^^l,,{it,it.4)- (188) 

For t(^^^) the jet functions are C(A"^). 

Finally, for T^'^^"'^ the time-ordered product has a structure that will produce the same soft 
operators, O^^ and P^f , as T^^-^^ and thus has shape functions ft\K^^2) fi'^\ii,^t)- 
The jet functions for T^^^") are 0{X-^). 



F. Summary of hard coefficients 

To determine hY'^^ we need to compute the traces in Eqs. ()144|1 . ()158j] and ()17()j) . Since 
there is no possibility of confusion, we shall now drop the bracketed factors, [2—], in the 
superscripts. For B X„ez/ {f = u), expanding to leading order in yn ^ uh gives 
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For B — i> Xs7 (/ = s) we expand around xh = 1, which gives 
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VII. SUMMARY OF DECAY RATES TO NLO 
A. Discussion of NLO results 

To facihtate the computation of the NLO corrections to the decay rates, they were divided 
into several pieces, as discussed at the beginning of Sec. IVll The full results at order A/mf, 
and all orders in as are quite complicated. However, now that factorization for the NLO 
rates has been achieved, we are able to expand consistently in factors of evaluated at 
perturbative scales. In this section, we discuss which terms are kept for phenomenological 
purposes. 

(2) f 

The most complicated contributions are the NLO correction to the Wi, namely , 
where / = s for 5 ^ X,7 and / = m for B XJu. From Eq. (fTTH|l in Sec. ED the NLO 
factorization theorem is 

'imf, Jq 

+ E 

r=l 
4 

r=3 
6 

r=5 ^ 

+ ... , (191) 

where j = 1,2 and j' = 1, 2, 3. The ellipses denote terms given in Eq. ()118|) that have jet 
functions that start at one-loop order. These additional terms are of order 

(192) 

71 mb 







h'/{n' 


■P) 


mb 
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where /xq — -^/mfeAQCD, so they are suppressed relative to the terms displayed in Eq. (jiyij] . 
They also induce dependence on new unknown shape functions, go-26, which makes it pro- 
hibitive to include them phenomenologically. Thus, throughout this section we shall stick to 
the terms displayed in Eq. (jl9ip . which enter at order A/m^ (r = 0-4) and 47ras(/io)A/m6 
(r = 5, 6). For convenience we drop the subscript 1 on i//^^'"^^ in this section, since there is 
no possibility of confusion. At tree level, the jet functions were computed in Sec. I VI CI and 
are 

= 6{k+) , J^-^\kl) = ^^''ll -_ , (193) 

The result for JT'-'^^ at one- loop order can be found in Eq. (j93|l . Examples of diagrams that 
contribute to Eq. (jl91|) can be found in Table IIVI Results for the hard functions h^/ were 
given m Eqs. (fTUni) and (fTTIT| and for h]^'^^ in Eqs. ffM and (IMUl . The running of the 
/i[^'s in Eq. ffT!ITl) is the same as the running of the /i-^'s in Eq. (HUO)), and the /i.-^'s depend 
only on the running of the Cj's computed in Ref. 0|. The shape functions /q^q are defined 
by Eqs. (fTT^ and (fTT^ . 

More 0{\^) power corrections occur when one switches to hadronic variables in the LO 
factorization theorem; these were given in Eq. pi3|) . For phenomenological purposes terms 
that start suppressed by the factor in Eq. ()192j) are not included, so we include the (A1 + 3A2) 
term but not the A term. 

Finally, the simplest source of power corrections is the kinematic expansion of the pref- 
actors in the decay-rate formulae, which depend on which rate we consider (see Eqs. ()42|1 - 
(j44p ) . In addition there are corrections from the expansion of n-q in the /ij's; these give 
the /if^'s in Eqs. ffTTa and dTUHl) . Both of these sources involve the same jet function and 
shape function as the LO contributions but different hard coefficients, which we shall denote 
by , and G*-^ in the next section. 



+ 



Mm) (ktm) (ktm: 



B. NLO Results for the Endpoint Decay Rates 

In this section we combine the pieces to arrive at the NLO results for decay rates. The 
hadronic dimensionless variables that we use are 

Xh= — -, xh= , yH=^, sh= yH= , uh= , (194) 

niB ttlb ttlb ttlb 

where only the first one is for B X^'-y and the rest are for B X„£z/. Note that 
sh = uhVh uh = {^~UH){l — yH)^ so these variables are not independent. The decay 
rates derived with SCET are valid for spectra dominated by the endpoint or SCET region. 



If one desires, he or she can straightforwardly treat xh in the doubly and triply differential rates in a 
different manner from what we have done here, by using Eqs. (|34|1 and (|35|l instead of Eqs. 144|l and (|43|l 
and making the desired expansions. 
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in which uh/Uh ~ -^1/ 1 ~ ^li ~ ^^h-> where \\ is a small expansion parameter which 
one can choose to be ~ 0.2 or ~ A/m;, ~ 0.1. In this region we write the decay rates as 

idv /irfr\^^ /idr\'^^° , , 

+ +■■- (195) 



VdZ \VdZJ \VdZJ 

where Z denotes a generic choice of phase-space variables. The doubly differential and 
triply differential decay rates that we consider were discussed in Sees. IIIBI and III (.1 We 
caution that if the SCET expansion is used for decay rates integrated over a larger region 
of phase space than the SCET region then it is important to check that sufficient smearing 
has occured so that the contributions outside the SCET region do not cause the power 
expansion to break down. In Eq. 1)411) we proposed a method one could use to combine 
results from the local OPE with those from SCET even if a larger region of phase space was 
desired (which become relevant when radiative corrections are included). Note that the \h 
parameter provides a means of testing for cases where a pure SCET expansion is valid. Here 
we present results for arbitary Xh-, and leave the investigation of the expansion in different 
decay rates to future work. 

For convenience, we define a shorthand notation for the convolutions that appear at LO 
and NLO: 

[J ® fWx.Px) = jdkl J{p^kl,fi) f{kl + A-p+,fx). (196) 

At LO we simply combine the hf'^^s and /i°"'s from Eqs. (jlOOp and (jlOlj) and find 

1 dV 



1 c/3r« 



LO 



Fg dx^* dyjj duH , 



F^ dyHduH 



HS [j(o) ^ j{0)] ^rnB{l-xl),mB) , (197) 
H^'iyn) [J^'^^f^%rnBUH,mByH), 



where the convolution of the LO jet function and shape function is given by Eq. ()196|) (with 
one integration variable k~^) and the hard coefficients are 

H'i^^) = mB[Ci'^-lci'^-Ci'^]\ (198) 
H^{fx,yH) = l2mB yHa-yn) [c[^^ ]\ 

Note that we have suppressed the functional dependence of the coefficients, which once 
perturbative corrections are included is Cf\^,mB) and Cl^\fi,mB,yH)- Our W^^^^s for 
both B — s> Xuii^ and B —* Xg'y, and hence the corresponding LO triply and singly differential 
rates, respectively, agree wth Ref. 0. These results also match those of mM 
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At NLO, combining all contributions, as discussed in the Sec. lVII leads to the following 
expressions: 

+ (l-x],)G« [j'°'®/<»>](mB(l-iJ),mB) 

6 

+ Y^H^i^B) [:r("^)®/,](mB(l-x^),mB), (199) 

j=0 



6 

+ J2 Hjiyn) [J^"^^ ® /,] {rriBUH, mByu) , (200) 

j=0 



+ unG^'iyH) [j(°^®/(°T(mBMH,mByH) 

6 

+ 5^ Hfiyn) [J^""^^ ® /,] (msMH, mByu) , (201) 

j=0 

where hq = rii = n2 = 0, = = —2, and = uq = —4, and the notation [jT'*-"^-' ® /j] 
denotes the convolutions displayed in Eq. ()191|) (see also Eq. ()114|) ). From the NLO phase- 
space factors in Eqs. (j42p . (j44j) and (j43|) . and from the NLO terms in the expansion of the 
hi, i.e. Eqs. XTU^ and (fTTI^ . we find that 

G'{m,) = -SmsicP-^ci'^-C^y, (202) 
GUvh) = 12ms{yJ, - l){Ci"^)\ 

GUvb) = 12m^{2(l-y^)(C;'^))^-[C7;^)+Cr^ + ^CrT}, 



G^{y^) = 2mByH[(2y^-6)(Ci''^)^-yH(y^ + 3)(C;'^) + C('^))d'') 
-yliCi'^r - 2(l+y^)(2Cr + C7r)Cr 
By using Eqs. ()189|1 and ()190|1 . we find that the hard coefficients for the terms with subleading 
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shape functions are 



TjS TjS TjS 1 TjS 1 TtS 1 Tja r) TTO fl 



Ho = HF = tni^ - 2^h), HF = 2yH{2yl + - 2), //I" = - 6^^, , 

HF = 2yH-Ayl, HF = 8yl-8yH, = Ay^ . (203) 

The Hi_Q factors displayed here will be corrected by as{mi,)/TT terms at one- loop order. 
Note that the shape functions /s^g depend on the flavor of the light quark in the four-quark 
operator and therefore differ for B ^ Xg and B — >■ X^- 

Equation ()199p describes B — > in the endpoint region. In Eqs. ()200p and ()20ip . 
convergence of the SCET expansion requires uh ^ Vh- For Eq. ()200p . we have in addition 
made a cut on xh such that 1 — xh ~ A^. In contrast, for Eq. pOlll the full range of xh 
has been integrated over. We can also straightforwardly obtain the rate (PV^ / dq^dm\, or 
equivalently d^V^ / dsndyn, by changing variables from {yfj^uu} to {yn^ sh} in the Wi and 
re-expanding. This is done by using yu = C — sh/C ~ ■ ■ ■ and uh = sh/C + ■ ■ where 
( = 1 — yn + sh and the ". . ." terms are not needed at NLO in the power expansion. 
The result may then be substituted into Eq. ()45|1 . Integration over yn then gives the 
spectrum, dT^/dsn, for which restriction to the endpoint region requires cut{sH,yH) = 

{sH/e<xi/{i+xir}. 



C. Singly Differential Spectra 



In this section we give results for singly differential spectra at NLO that are useful for 
measurements of \Vub\- Since our goal is to explore the effect of the sub leading shape func- 
tions, in this section we shall work at lowest order in ag for the hard and jet functions, both 
at LO and NLO. For the shape functions we shall use f[^\uj) = 2ujf^^^{uj) jsj] and define 



(204) 



where a prime denotes a derivative, as well as 



6{kt) - 5{ki 



dk-^ d}\i<2 dkr^ 



mm) (ktm) (Km 



-n'6{kt)6{k+)6{kt 



x/(S(A;+ + A-p^ 



(205) 
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We shall put an additional superscript s or m on F^^Q^p'^) to distinguish the original four- 
quark operator with strange- and up-type quarks. For B Xs'j the rate dT^/dxJj in the 
endpoint region is equivalent to making a cut on xjj. The necessary results already appear 
in Eqs. fll97p and p99|) and combining them gives 

1 dV^ — 

--— = mB[l-3{l-xl)]F{mB{l-xl)) + [mB{l-xl)-A]F{mB{l-xl)) 

+F2{mB{l-xl)) -F3{mB{l-xl))+F^{mB{l~xl)) 

-8naM F^{mB{l-xl)), (206) 

where 1 — ~ A^. Here the —3(1 — xjj) term comes from expanding the kinematic 

— (2) 

prefactors, while the term that depends on [ms(l — xh) — A] comes from reducing fl to 
/(°) (and we write it in terms of F for convenience, even though this includes some formally 
higher-order terras). The rate in Eq. (j206p for B — > Xg'j has been computed previously 
at NLO in Ref. [ij. The relation of our shape functions to those in Ref. js^ is shown 
in Table IVIII Our result in Eq. ()206|) agrees with theirs, up to the F^^^ term, which was 
computed here for the first time.^^ 

To compute dV^/dxH for B — > X^/z/, we can use Eq. and case i) of Table HTl Treating 
= 0.026 ~ A^, we can set = even at NLO. To ensure that the SCET expansion 
converges requires that one makes a cut on uh- First consider the xh spectrum. Since 
uh < ^ — Xh, making a cut xh > xjj restricts uh to the desired small values. However, at 
NLO accuracy this is not equivalent to our uh < ^^hVh definition of the SCET region of 
phase space, and an additional term, K{Xh,xh), is added to correct for this. Depending 
on the choice of other cuts, the error in including a larger region of phase space may be 
power suppressed. The parameter Xh provides us with a way of testing this by comparing 
Xh = 0.2 and Xh = 1- We present our final results in a manner that makes it easy to take 
the Ah ^ 1 limit for situations where a large enough region has been smeared over that this 
is the case. For the xh spectrum the result is 



1 c/F" 



F^ dXH 



VH 



2 j duH 1^5(1 - AuH)F{mBUH) + {A-niBUH) F^mBUn) 



- F2{mBUH) - SFsimBUn) + ?>Fi{mBUH) - 247ras(/^o) F^imBUn) 
+ KiXH,XH), (207) 



In making this comparison, note that Ref. (51| gave their result with partonic variables, whereas we have 
expressed ours in terms of hadronic variables here. 
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where 



K{Xh,xh) = 2 J duH {2uh-SXh) mBF{mBUH) + {A-mBUH)F{mBUH) 

{1-xh)XI 

-F2{mBUH) - {uh-2\b) 3F3{mBUH) - SF^^itibUh) + 2'iTcas{iJ^o) F^imBUn) 



niBFimBUH) \. (208) 



The A//-dependent term, K{Xb,xh), arises when one writes the integral over yu: as 
Ji_^^dyH + dyu^^ ^fr- — 1 + Xh^, which is done to ensure that Uh < ^'hVu- 

ter inputing models for the shape functions one can use the K{Xh,xh) term to check the 
consistency of the operator product expansion with this level of smearing. 

The rate dV^/dxn has been computed at NLO in Ref. [s^ and we can compare our result 
with theirs. To do this we take —>■ 1, which sets K{Xh,xh) —>■ 0, since the restrictions 
we imposed on the phase space were not considered there. We also must convert back to 
partonic variables, which means dropping the (A1+3A2) term in Eq. ()204|) . After doing this, 
we find agreement with their Eq. (35) on the coefficients of the LO F term and the F2 
and F4 terms, but we disagree on the uhF term and the F3 term. (Again the F5 term is 
computed for the first time here, so no cross-check on this is possible.) The coefficients of 
our F3 term also disagrees with the RPI constraint derived in Ref. |83], which predicted 
that they occur in the combination tubF — F^}^ We found that this combination occurs for 
B —>■ Xg'-f, but not for B X^iu. In the next section, we present a non-trivial cross-check 
of the coefficients of our result in Eq. ()207|) . namely that when re-expanded it correctly 
reproduces terms in the local OPE up to 

Next we consider the Px spectrum, namely dT/dyjj. Integrating the doubly differential 
rate from Eqs. p97j) and ()2()1|1 with a cut on uh, i-e. uh < ^jjVH, restricts one to the 
triangular SCET region shown in Fig. |21and gives 



1 rfE" 



dyn 



= 2 du 



H |mB[(3y|-2y|^) +2uH{yH-SyH)]F{mBUH) 

+ (A-mBUH){yH-2fH)F{mBUH) + {2y]j + yH-'^yH)F2{mBUH) 
+ {2yjj-3yH)F3{mBUH) + {yH-2y]j)Fi{mBUH) 

+16TTas{fJ^o){yH-yH)F^{mBUH) + 8TTas{fio)yHFQ{mBUH) } ■ (209) 



Another possibility is to consider the spectrum, which is dT/dun, where we now integrate 



We found that they have an overall 2 typographical error in this equation (see also j56l]'). 

Note that, from the point of view of the factorization theorem, such a reparameterization constraint would 

be very interesting, since it would give a relation between the jet functions j'^^'^^ and j/'"' to all orders in 

OLs, even though these operators appear to be defined by unrelated matrix elements, Eqs. H83|l and l|152(l . 



61 



with the cut on yu-, obtaining 
1 rfr" 



To duH 



\mB [{l~2u^+u^) - uuiJi-Qu^+u^)] FimBUH) 



(210) 



o o 

-i(5-9M2+4M3) F^{mBUH) - \{l + ?>u^ -Au"") F^imBUn) 
o o 



[l-3u' + 2e) F^iruBUH) + 87ra,(/io) Fe"!^^^//) 

d 

where m = is a parameter of order 1 for A^/ ~ 0.2. Given sufficient smearing, one 

can take — > 1 and the u terms become subleading. The subleading terms in dT^/duH 
provide power corrections to the phenomenological analysis using this spectrum in Ref . [Uj . 
We leave the derivation of dV^JdsH with phase-space cuts to future work, and so have not 
compared this rate with Ref. |5q . 



D. Comparison with the local OPE 

A non-trivial check on our power-correction results can be obtained by expanding the 
subleading shape functions in a manner appropriate to the case where the local OPE is 
valid. Essentially this means expanding 

5{t - m-D) = 6{e+) - 5\t)in-D + ... . (211) 

The comparison will be made at the level of the partonic dV^/dx decay rate, and we define 
w = mBUH — A. From Eq. ()207p we have 



f ^ dx 



2(—\ [ "duH {niBil - 4uj^)/(°)(A-mBU^) + i/(2)(A_^^^^) (212) 

- (^i + 3^^) /(oy(A-^^^^) + (K-rriBUH) f'\'K-mBUH) - fi'\A-mBUH) 

- SFs^uibUh) + 3Fi{mBUH) - 247ra^(/io) F^{mBUH)^ + K{\h,xh) 

2-^ / dw\l + /(°)(-u;) 

— /f (-^) - —Fsiw + A) + —F,iw+A) - ^^I^lM F-iw+A)] + K{\h.xh), 
ruh nih ruh nih J 



where x = 2Ei/mb and f g = rn\/m\ Fq. In writing the second equality, the (Ai -|- 3A2) term 
was cancelled by the change in the upper limit of integration. The — 4A/mb term cancels the 
leading term in the expansion of {mB/mb}'^, while the terms that cancel higher terms are 
beyond order A^ in the SCET result. The F5 term gives an as/mf term, which we shall drop 
below. In expanding the shape functions, we obtain singular functions peaking at x = 1, so 
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it is safe to take Aw = 1 and drop K{Xh,xh)- For the remaining shape functions, the local 
expansion gives 

f'\-w) = 6{w) - ^ 6"{w) - 6"\w) + ..., 

o 18 

/o^'^(-^) = -(Ai + 3A2)5'H + ^5» + ... , 
Fs(A+w) = -^Ai(5'H + ... , 

F^iA+w) = -A2 6'{w) + . . . , (213) 

where Ai, A2, pi, P2 are the standard HQET parameters, which are matrix elements of lo- 
cal operators, r is a combination of time-ordered products of HQET operators, and terms 
beyond l/mf have been dropped. Substituting this result into our Eq. ()212|) and integrat- 
ing by parts using the relations js^] / dw[w6"{w)] = J dw[—26'{w)] and J dw[w6"'{w)] = 
J dw[—36"{w)], which are valid with smooth test functions, we find 

+^ 6'{1 -x)--^ 6"{1 -x)-^ 6'{1 - x) - 6'{1 - x) . 
2m'l 9mi 3m^ 

This agrees exactly with the result obtained from the local OPE in Refs. jisi Eo^ . (Note 
that we do not compare the 1/ml annihilation term, which does not arise from one of the 
shape functions appearing at order A^ in SCET.) 

Of the terms in Eq. ()214|1 . it is those proportional to Ai and pi that test the difference 
between our results and those in Ref. [H^. Ref. [H^ also obtained the Ai result in Eq. ()214j) . 
and even though we disagree on the unf^^^ and F3 terms, the combination of the two gives 
the same Ai result. For the pi terms, their pi5" term agrees with Eq. ()214p . but the pi5' 
term does not. In the very recent paper [3] it was pointed out that from the local OPE the 
coefficient of the pi5' term should be —5/3, as in Eq. ()214j) . rather than the —1/3 quoted 
inS. 

Ref. [8J| went further to advocate a different approach to the shape-function region that 
involves using an unexpanded 6-quark field, and doing this obtained a result with subleading 
shape functions whose expansion is consistent with the local OPE in Eq. ()214|1 . However, 
their result for the power-suppressed dT^/dx decay rate does not agree with the rate obtained 
here. In particular, our —4uhF term is not present there, and instead of our — 3-F3 term 
they have the result "4Fi + 2G2 — 3G3". Their "G2" term is defined by operators that, in 
our analysis, can only show up suppressed by at least one factor of through jet functions. 
In fact, the operator structure of our result actually agrees with the original one in Ref. js^l, 
rather than the one in Ref. No proof of factorization has yet been achieved for this 

approach with the unexpanded 6-quark field, and it is conceivable that this may help to 
explain why a different structure of operators was obtained. 
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VIII. CONCLUSIONS AND DISCUSSION 



In this paper, we have computed a factorization theorem for the leading-order power 
corrections to inclusive B — »• Xgj and B — »• X„£z/ decays in the endpoint region, where the 
X is jet-like. In particular, we have shown that these power corrections can be fully cate- 
gorized and thus treated in a systematic fashion using the Soft-CoUinear Effective Theory. 
A main result of our analysis is that perturbative power corrections to the decay rates can 
be systematically computed, and our result explicitly disentangles hard factors of as{ml), 
coUinear jet-induced factors of as(m|^), and soft ("q;s(A^)") non-perturbative QCD effects. 

In addition, our results can be used as a starting point for the systematic resummation of 
Sudakov double logarithms in the power corrections. To achieve this, one needs to compute 
the anomalous dimensions of all the operators we have defined that appear in the subleading 
factorization theorem. Some of the terms here are already known. In the body of the paper, 
we have shown that if we consider only subleading terms with non-vanishing jet functions 
at lowest order in as, then the logarithms that can be resummed into the hard function in 
these NLO contributions are identical to the analogous logarithms in the LO result. (These 
logs can be thought of as occurring between the scales ml and m^A.) There are additional 
logarithms that are sensitive to the split between the jet and soft functions (logs between 
m^A and A^), which require knowledge of the anomalous dimensions of subleading soft 
operators. The latter are very unlikely to be universal, and have not been computed here. 

Our main final decay-rate formulae have been collected in Sec. IVIII At lowest order in 
as, they include a derivation of the power corrections for the triply differential B X^ii^ 
rate. Results have been derived in the literature for dT/dE^ in B Xsj jHH and the 
singly differential B X^ii^ rate dT/dEi 5^ 5^, 8^ (and dV /dm\ jH^), and a comparison 
was given in Sec. IVII CI Agreement was found for B — » -^s7; but for dT/dEi we found 
disagreement on two terms at subleading order. (A check on our dT/dEi result was obtained 
by expanding it to compare with the local OPE, and full agreement was found up to l/ml, 
as discussed in Sec. IVIIDH Using our results, we found it straightforward to present the 
power corrections to doubly differential rates, as well as other singly differential rates such 
as dT/dp]^ and dT/dp~^. 

On the phenomenological side, our most important result is the identification of two new 
shape functions, which involve four-quark operators and have not been previously considered 
in the literature. They are denoted by f^^ (-^s.e), and definitions can be found in Eq. ()133|1 . 
In the endpoint decay rates they induce power-suppressed terms, which are quite large, of 
order 

A , , 

47ias — . (215) 
nib 

Since Anag ~ 4, these power corrections might numerically dominate over those that are 
simply of order a^A/rrii,. We have given results for the effect of these shape functions in 
all the considered decay rates. In our results for the decay rates, the numerical prefactors 
for /5,6 turned out to be sizeable (e.g. 24:7ias for /s in dT/dEi), which justifies including 
the factor of 4 in Eq. (|215p . For the extraction of \Vub\ from dT/dEi, the important thing 
to consider is the difference between how these new shape functions affect B Xud-v and 
B Xs7. In this case, comparing the combinations of F and F^ in Eqs. ()2U6|) and ()2U7p . 
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we see that the mismatch is 



(216) 



where the index m or s denotes the fact that these shape functions involve different flavors 
of light quark. To obtain a numerical estimate we approximate / {rriBF) ~ F^ / {msF) ~ 
A/rrih ~ O.le' and find that they can cause a deviation of 



where e' denotes any additional dynamical suppression from the non-perturbative functions. 
Even for e' ~ 0.1 — 0.3 these terms provide a sizeable new uncertainty for the cut spectrum 
dT/dEi approach to measuring \Vub\- From Eq. ()210|) we see that the situation is only 
slightly better for a cut dV /dp\ spectrum. Our conclusion is that, until theoretical methods 
are developed to control these new subleading shape-function effects or better bounds on 
these contributions are obtained, it will be challenging to argue that methods sensitive to 
the shape-function region are trustworthy for inclusive determinations of \Vub\ at the ^ 10% 
level. One possible future direction is to derive experimental bounds on the new subleading 
shape-function effects by comparing endpoint-dependent methods with different spectra and 
different cuts on the phase space. It would also be useful to find model-independent ways of 
determining the size of the subleading shape functions that go beyond the simple dimensional 
analysis used here. 

Theoretically, there are several avenues for future work on B ^ and B —>■ X^ii^- 
These include the calculation of perturbative corrections in the factorization theorems at 
subleading order, as well as a complete resummation of Sudakov logarithms. It would also 
be interesting to consider the structure of the subleading factorization theorems in moment 
space, as opposed to the momentum-space version considered here. Starting with our triply 
differential B Xu^v result, one could derive other doubly and singly differential decay 
spectra and consider their phenomenological implications. More formally, it remains to be 
checked that the convolutions that appear in our subleading factorization theorems actually 
converge when the functional forms of the jet functions are considered at higher orders in 
Us- From a formal standpoint this is necessary for a complete "proof" of these results as 
factorization theorems. However, from a pragmatic standpoint this can be checked as each 
new phenomenologically relevant term is computed. We are not aware of any factorization 
formulae where convergence problems occur at higher orders in the perturbative expansion 
of the kernels when they are not present in the leading non-vanishing kernel results (the 
convergence of which we have checked). Finally, it should also be possible to extend the 
techniques used here to closely related physical cases such as deep inelastic scattering for 
X — > 1 (i.e. Bjorken a; ~ 1 — A/Q). 

Note added: Recently the paper appeared, in which a different type of power- 
suppressed correction was studied near the endpoint region, namely those corrections sup- 
pressed by A/ A, where A depends on the cut on the photon energy spectrum. These 
corrections involve using the leading-order shape function and give a measure of the effect 
of extrapolating away from the SCET expansion region we discussed. The results derived 
here should enable this analysis to be further extended to include the power corrections 
suppressed by a factor of K/mb- 



167ra,(/io) O.le' - (180%)e' 



(217) 
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APPENDIX A: EXPANSION OF THE HEAVY QUARK FIELD AND DERIVA- 
TION OF THE CURRENTS 



To derive the power expansion of the heavy-to-light current J = ip^'^^T^^'^^ at tree level 
we just need the expansion of the light- and heavy-quark full-theory fields, ip^'^^ and tp^''\ in 
terms of SCET coUinear fields and soft heavy fields (h^) respectively, ip^'^^ also contains 
terms with usoft light quark fields, but they do not contribute for inclusive processes until 
higher order, and so we neglect these terms here. For the light quark field we have 0| 

— AW-^ + WiB^ W^)^ 



(9) 



1 



in -De 
1 



in 



(Al) 



■p2 



Here we have used the collinear fields ||75| 

in-D = in-d+gn-Ac+gn-Aus , iD^ = iD^ + WiD:l;^W^ , 

in-D = in-Dc+Win-DusW\ (A2) 

which obviously will give an expansion with terms that are individually usoft and collinear 
gauge covariant at each order in A. 

For the heavy quark fields we start with the QCD Lagrangian C = tjj^''\ip— nib)^^''^- We 
then divide the quark fields into on-shell and off-shell terms, ■ip^'^^ = exp{—imv-x){hv +V'-b), 
where -^h^ = h^. Unlike in HQET, the off-shell field ipB does not satisfy ■^tpB = —^^b, since 
collinear gluons give off-shell quark fields that have quark components, ■^ipB = '4'b- This 
gives the Lagrangian 

C = KivDusK + i^B[ip+ m^^i - \)\i\)B + K{ip- Wb + ^b{%P- f)K , (A3) 

where the subtraction in (ip—'p) removes the terms that vanish by momentum conservation. 
Varying with respect to ipB gives 

m,(^ - 1)]^B = -{ip- f)K , (A4) 

which is a higher-order version of the LO equation for ipB derived in Ref. |37|. To solve 
for 'i/'B in terms of at higher orders we use a strategy proposed in Ref. [73[, namely 



expanding ipB = "^Pb^ + "^Pb' + ■ ■ ■, and considering the solutions order by order in A (noting 
that ~ A^). To facilitate this we rewrite Eq. ()A4|) as 



.,(4) 



A(o) 



(3) 



+ ...) 



gn-A, + g4j + A^^^ 



(A5) 
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The complete set of equations to solve is then 



-^gn-AcK 



(A6) 



(n) 



{2),/,(n-2) 



where n > 6. Obviously, the crucial point is to invert the operator A'^*'^ When acting on a 
collinear field, the inverse [A*^°)]~^ is simply given by 



[A 



(0)1-1 



[l + f) 
2mhn-v ' n-vin-Dc 



+ 



(A7) 



Owing to the l/{in-Dc), this solution is ill-defined when acting on an ultrasoft field and will 
not suffice to solve for ip^^-^K For the solution to the first equation we find 



in-D, 



(A8) 



which, as expected, is in agreement with the LO solution in the appendix of Ref. |37|. Thus, 
at lowest order the solution for the full-theory field is ip = hy + ip^ = Wh^. At the next 
order, we find that 

(A9) 



1 ^ 

rubn-v 2 



W 



■^pW^igp'^W 



h„ - 



-W 



n-v 



-^WkgvBlW 



where 



igB^j^ = [in-D'',iD^^']. (AlO) 

The result agrees with Refs. 7^ 2^, except that we have written the last term as v-B^^ to 
emphasize that it must start with at least one collinear gluon. 
To proceed we let 



and then write the remaining equations for ijj^^^ . Using the equation of motion for W gives 



(All) 



A(o) = |p+m,(|^-l), [AW];1 



2mi,n-v n-vV 

and so, after we multiply on the left by W\ the remaining equations become 

A(2) - ^' 



(A12) 



(A13) 



'-{W^in-DW -in-Dus) 
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where n > 5. 

To solve Eqs. (fAT8jl we invert A'-^''. For terms with at least one collinear gluon field, 
such as A^^''^\ we can use the inverse [A^^^]"^ from Eq. ()A12|) . Only terms that are purely 
usoft, including ipushv and possibly ipus'ijjB ~'^\ need to be handled with care. For these 
terms a purely usoft ip^^'^ suffices, in which case A*^'') = m{'^ — 1). The subset of purely 



usoft terms in Eqs. ()A13j) is given uniquely by the terms m(|^ — 



g) 



-ipushv and m{j} 



= —ipusi^^B which are the same as we would find in pure HQET. Thus, the 



solution for the pure usoft terms is 

1 



us 
B 



2mb + iv-Dii 



1 {i-f) 



2mh 



(A14) 



where it is sufficient to use this [A'^^^J^j^ since (1 - i))/2%l)l' = ifj^'. From Eq. ^KT^ . the 
solution for ip^^^ is then 



(A15) 



W- 



2m, 



Arribn-v 



■:pW^ign-BM 



hy-W 



n-v 
n-vV 



^WHgn-B,W 



-W 

+w 
+w 
+w- 



rribn-v iV -I LP 



nib {n-vY 



-^WUvD^W 
V 



hr. 



rrib {n-vY 



^WHgp^w] l^WUgvBlW 



in-v] 



V 



-^W\gvBlW 



From the expansion of the fields we obtain the expansion of the currents at tree level. 



^(3) 



(4) 



(A16) 



The leading-order SCET heavy-to-light current is 
For the currents suppressed by A we obtain 



(A17) 



in- Dc 



-rwL, - 



n-vm ' 



2in-D, 



-igp[WK 



;igvBlWK. 



n-v (in-Dc)'^ 

Making the field redefinition and putting in the most general tu, dependence consistent with 



68 



RPI gives VI 



17 



= l{i.i^iW)^^T{Y^K), (A18) 

where B-^ is defined in Eq. For the currents suppressed by we find 

J(2) = lj:w'^K - —InV (^ + -^zgn-MWK (A19) 

2m n-v \4m in-DcJ in-Dc 



[n-v)'' tn-Dc [in-Dc) 

[n-vym \ m-Dc m-D^ J m-Dc 



where 



ign-M = [in-D^, in-D]. (A20) 



1^ Note that for j'-^"^ we find a field strength B^^ rather than a v-D']_ acting on W. 
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It will be convenient for us to label the terms as follows: 



j{2c) ^ -"^Ij: _ \ ign-MWK, 
n-v \in-Dc) 

J^''^ = —^nr'^-^^9n■MWK, 

n-vm ^in-Dc ^zn-iJc 

j(2/) = ^_^^Y-^tP:,-^tg^_,WK, 

n-vm tn-Dc tn-Uc 

J^''^ = . \, ^nT-^zv-Dji-^zgPiWK, 
[n-vym tn-Dc tn-Dc 

[n-vym in-Dc [m-Dc) 

J^''^ = -^^nT-^tv-Dj _ ^ igv-BlWK. (A21) 
The terms J^^a) j(2/) 

agree with those found in Refs. BEi, 

which use the position- 
space formulation of SCET. J*^^"\ J*^^^) and J'^'^^^ correspond to the third, last and fourth 
terms, respectively, of Eq. (41) in Ref. 0, while J*^^'^-', J^'^^'> and J^'^^'^ correspond to the 
second, last and third terms, respectively, of their Eq. (43). The rest of their second-order 
terms correspond to Lagrangian insertions in our SCET calculation. The terms J^^a) to J^^i) 
do not appear in Refs. [g^ ll^, because they set f ^ ~ A (more generally, the assignment 
f ^ ~ 1 used here is allowed). We follow the common practice of dropping the v^ terms in 
our analysis of the factorization theorems by picking a frame where = 0. 



APPENDIX B: REPARAMETERIZATION INVARIANCE FOR THE CUR- 
RENTS 

The choice of the light-cone vectors n and fi must satisfy ra^ = = 0, n-n = 2, 
but otherwise is arbitrary. Therefore one can make the following reparameterization (RPI) 
transformations 

(I) / ^ + ^jj^ j ^ n^ ^jjj^ ( n^ ^ {1 + a) n^ (Bl) 



^n^^n^ in^^n^ + e^ \n^ {I - a) n^, 

where A-*- ~ A, while ~ a ~ A*'. 

It is useful to separate the RPI transformations according to the power suppression in A 
that they cause. We shall denote this by a bracketed superscript. For type-I transformations 
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we have, at the same order in A, 

6i{n-D) = 6j{W) = 0^ 



A 



>in.D) = A-.Dt , 5r '(D-) = -^n-D, 



M 

c ; 



Sf\^n) = 0. (B2) 



Type-I transformations one order higher in A are 



r'te.)=4^. (B3) 



For type-II transformations we have, at the same order in A, 

Snin-D) = 0, 



= Sii"\w)=0. (B4) 
Type-II transformations one order higher in A are 



n- D r ^ 



n-Dr 



(B5) 



where the differential operators do not act outside the square brackets. Furthermore, 



= LWpiW^^^-^ . (B6) 

Note that Table I of Ref. js^ seems to suggest that the transformation of W involves the full 
covariant derivative with both coUinear and usoft parts. However, the Wilson line referred 
to in that table is constructed from n\A(. + A^^] for the Wilson line built out of n-A^ alone, 
the equation of motion is in-DcW = 0, and hence Sii{W) = 6ij\w). Nevertheless, there is 
a potential complication related to the choice of the full covariant derivative. One possible 
reparameterization-invariant choice is in-D = in-Dc + in-Dy_s, iD^ = iD^ + iD:^^. With this 
choice, one must redefine n-Ac and A^ if one wants the Lagrangian to be invariant order by 
order in A under coUinear gauge transformations (75|]. This turns out to be equivalent to the 
choice given by Eq. ()A2|) . with the covariant derivatives defined in terms of the new gauge 
field. Furthermore, the Wilson line built out of the new gauge field has the same equation 
of motion as previously, and hence the same transformation law. 

Using the above transformations and the definitions ()A10|) and ()A20|) . we obtain 

6i''\zgBl^) = -^n^Ai[zn-D,,zDi^] + [6p{WzDi^W^),tn-D4 

= -^mgA^-Bl, (B7) 
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since the second commutator in the first hne above equals [WS[^ \iD:^^)W^ ,WVW^ 
W[5i^°\iD^^'),P]W^ = 0. Similarly, we find that 



sf\igBf) = , Sf\ign-M) = igA^-B^ . 
Then, it is straightforward to show that 

It follows that J(i) + J(2) = 0. 

Repeating the process for type-II transformations, we obtain 



(B8) 



(B9) 



n 



= ^gei.G^';-^lgn■M, 



(BIO) 



where gG^l = i[iD^i' ,iD^% Similarly, we find that 



5[f\igBf) = '-ige^-Bl , ^^^(i^n-M) = . 



(Bll) 



Then, 



in-Dn ^ 



2 in- D, 



2 in- D, 



+ ^^^^j^^^is^.D,^^^WK + ^^^^i^^^^ 



1 



in-D, 



ie^-D^W 



2 ;r ^ 1 

+ — Cnr- 



n-v^" (in-Dc)^ 



ie^-D^- — —igvBlWK - v^e^igC^^WK 



+ igv-B1 



——ie^-D^W 
m-Dr 



in-Dc 

2 - 1 
hn, 1 irJ^ ._ „ iv-Dt 



n-v in-Dr 



{in-Dc) 



jige^-BlWK 



-^ni^j ^^^T ige^ - BIWK - 

^ in-Dc i'^^-^cr 

Now, since {in- Dc)~^ige-^ - B']Why — e-^-{iD^W — 'P±)hy, we can write 

1 



(B12) 



ie^-D^W 
in-Dc 



K, = -^e^-{iD^W - Va.)K = .J ge^-BIWK - (B13) 



in-Dr 



{in-Dcf 
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Substituting this into Eq. ()B12|) and combining the first and third terms, we obtain 

+— ar .., (te^-D^v-itD^W - V±)K 
n-v [in-Dc)^ \ 

-{ie^-D^v{iD^W -Va_)- {ivD^ie^-D^W - ie^-D^vV±)}K 
-ivD^e^-iiD^W -Va.)K 



in-Dr J ^ in-Dr 



c 

2 , „ 1 



■ <nr—^itv-D^e^-V±-te^-D^vV^)K 
n-v [in-Dc) 

= 0. (B14) 

The first term of the second-last equahty above is zero, owing to the equation of motion 
obtained from Eq. (j55j) . In the second term, the expression in brackets begins at one collinear 
gluon, ensuring that the 1 / [in-Dc] does not cause any singularities. Thus the label operators 
make this term zero. 

Using completely analogous manipulations, we also find the closure relation that involves 
the J^^^^ current. Thus the type-II constraints are 

) J(0) + ^(A^) J(la) + ^(A^) J(lc) ^ ^(A0)^(26) ^ ^(X^) ji2c) ^ ^(X^) ji2.) ^ ^(^ ji2,) ^ q ^ 



It then follows that ^i^'^jW +S^"'J^^^ + 5|f = 0. Dropping the j(29)-(2i) currents, 
these relations agree with Ref. j73l |. 

From the results in Eqs. ()B9|) and ()B15|) we see that the currents split into two sets: 

J(0) + J(la) + Jilc) ^ ji2b) ^ ji2c) ^ Ji2^) ^ ji2j) ^ 

jilb) ^ ji2d) ^ ji2e) ^ ji2f) ^ ji2g) ^ j(2h) (g^g) 

When we add non-trivial Wilson coefficients to all the currents, the number of Ui parameters 
they depend on is restricted by RPI, i.e. in the set with 7*^°^ they have one parameter u 
and in the second set with J^^'^^ they have two parameters, uji^2- In the second set the 
combination of fields restricted to have momentum ui and uj2 is determined by the manner 
in which the terms in the RPI-transformed currents cancel. Making the Y field redefinition 
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in Eq. ()57|) and putting in the most general Ui dependence consistent with RPI, we find 

J''"\'^) = ^ (Lw)^t[yHpIk) , 

where 

igB^j_ = [in-Dc,iD'^j_] , ign-Bc = [in-Dc,in-d + gn-An] . (B18) 
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